Critical adsorption and critical Casimir forces for geometrically structured confinements 



M. Trondle, L. Harnau, and S. Dietrich 

Max-Planck-lnstitut fur Metallforschung, Heisenbergstr. 3, D-70569 Stuttgart, Germany, 
and Institutfiir Theoretische und Angewandte Physik, Universitat Stuttgart, Pfaffenwaldring 57, D-70569 Stuttgart, Germany 

(Dated: August 14, 2008) 

We study the behavior of fluids, confined by geometrically structured substrates, upon approaching a critical 
point at T = T c in their bulk phase diagram. As generic substrate structures periodic arrays of wedges and 
ridges are considered. Based on general renormalization group arguments we calculate, within mean field 
approximation, the universal scaling functions for order parameter profiles of a fluid close to a single structured 
substrate and discuss the decay of its spatial variation into the bulk. We compare the excess adsorption at 
corrugated substrates with the one at planar walls. The confinement of a critical fluid by two walls generates 
effective critical Casimir forces between them. We calculate corresponding universal scaling functions for the 
normal critical Casimir force between a flat and a geometrically structured substrate as well as the lateral critical 
Casimir force between two identically patterned substrates. 



I. INTRODUCTION 

The confinement of a fluid which is close to its critical 
point at T = T c induces remarkable deviations from its bulk 
behavior. Typically, the boundary of a system affects the lo- 
cal structural properties of condensed matter within a layer 
of thickness of the bulk correlation length, which diverges 
upon approaching criticality as £±(t — ► 0) = £r |*|~ v , where 
t = (T — T c )/T c , v is a standard bulk critical exponent, and 
§q are nonuniversal amplitudes above (+) or below (— ) T c . 
For binary liquid mixtures, the disordered phase at T > T c 
is the one in which its two components are mixed, whereas 
in the ordered phase at T < T c there is phase separation in 
two phases being rich in one or the other component. (At 
lower critical points the disordered phase occurs for T < T c 
and the ordered ones at T > T c .) Accordingly, near the criti- 
cal demixing point T c the scalar order parameter is a suit- 
ably defined concentration difference of the two species, be- 
longing to the Ising universality class. The generic prefer- 
ence of a surface for either one of the two species forming 
the binary fluid amounts to the presence of an effective sur- 
face field which affects the order parameter profile close to 
the wall such that it is nonzero even in the disordered (mixed) 
phase for T > T c . Upon approaching T c , these so-called criti- 
cal adsorption profiles, describing the concentration enhance- 
ment near the surface, become long ranged due to the diver- 
gence of the correlation length. 1 - 2 3 In the semi-infinite geom- 
etry the transitions from the phase in which only the region 
near a single surface is ordered to the one in which also the 
bulk is ordered are known as the so-called extraordinary or 
normal surface transitions. 4,5 In the slab geometry the confine- 
ment of critical fluctuations of the order parameter gives rise 
to critical Casimir forces attracting or repelling the confining 
walls, which can be of significant strength. 1 - 6 ! 7 ' 8 ! 9 ' 10 This is 
the thermodynamic analogue of the quantum-electrodynamic 
Casimir effect originating from the confinement of vacuum 
fluctuations. 1112 Critical adsorption has been investigated ex- 
perimentally for flat substrates (see, e.g., Refs. ll3lll4lll5lll6lll7l 
and references therein). Also, there is strong experimental ev- 
idence for the occurrence of critical Casimir forces , 18 ' 19 ' 20 ' 21 
and recently a direct measurement of the critical Casimir force 
between a wall and a colloid has been reported. 22 So far the 



theoretical investigations of critical adsorption and the criti- 
cal Casimir effect have been focused on topologically flat sur- 
faces (see, e.g., Refs.l2ll 3 [l |l 3 23 24) or curved surfaces and 
colloids (see, e.g., Refs. I25ll26ll2l28ll29ll30h . 

Nowadays, experimental techniques are available which al- 
low one to endow a solid surface with a precisely defined ge- 
ometrical or chemical structure in the micro- and nanometer 
range (see, e.g., Refs. 31.32,33.34). Chemically and geomet- 
rically structured substrates are of major importance for the 
construction of a Tab on a chip'. 35 Critical adsorption and 
critical Casimir forces for chemically structured but topolog- 
ically flat confinements have been studied theoretically ! 36 - 37 
Typical man-made geometrical structures are periodic in one 
direction and consist of grooves with wedge-like shapes. The 
effects of such roughness and of the structuring of substrates, 
for example, on wetting phenomena (see, e.g., Refs. [38ll39ll40l 
and references therein), liquid crystal systems ^ 41 ' 42 ' 43 ' 44 or de- 
pletion interactions 45 have been found to be crucial. 

In view of this context, as a paradigmatic first step critical 
adsorption at a single wedge and a single ridge has been stud- 
ied at the critical point 46 and off-criticality;— It was found 
that, in agreement with earlier work^ the angle character- 
izing the wedge is a parameter which influences the critical 
behavior significantly. 

On the other hand, the effect of geometric structures of sur- 
faces on the quantum electrodynamic Casimir effect has been 
analyzed, too . 49 ' 50 ' 51 ' 52 ' 53 Periodic structures give rise to lat- 
eral quantum electrodynamic Casimir forces ^ 51 ' 54 ' 55 ' 56 ' 57 Re- 
cent studies have been focused on quantum electrodynamic 
systems in which lateral Casimir forces due to topological 
structures lead to a controlled motion. 58 - 59 ' 60 - 61 

The present work is supposed to extend these investiga- 
tions into various directions. As a paradigmatic model for 
a geometrically structured substrate, we consider a periodic 
array of wedges and ridges. One expects an interesting in- 
terplay between the externally endowed structures and the 
critical behavior of the system if the size scales of the sub- 
strate corrugation are comparable with the correlation length. 
Actually, this range is experimentally accessible because the 
correlation length can reach values up to several hundreds 
of nanometers (e.g., for a mixture of water and lutidine one 
finds £q" ~ 0.2 nm so that for t = 10~ 4 , corresponding to 
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T — T c ~ Q.Q3K in that system, one has t, ~ 70 nm). 22 ' 62 First, 
we describe the critical adsorption behavior of the order pa- 
rameter profile for the whole temperature range around T c in 
terms of universal scaling functions based on general renor- 
malization group arguments. We calculate these correspond- 
ing universal scaling functions up to lowest order, i.e., within 
mean field theory, and subsequently we discuss the experi- 
mentally relevant excess adsorption. Second, we study the 
critical Casimir effect for geometrically structured confine- 
ments and focus on the universal features of the normal and 
lateral forces emerging in binary liquid mixtures close to crit- 
icality. We calculate universal scaling functions for the critical 
Casimir forces within mean field approximation for identical 
chemical boundary conditions on both walls ((+,+) configu- 
ration). 

The paper is organized as follows. Section [TT] is devoted 
to the study of critical adsorption on a single geometrically 
structured substrate. General scaling properties of the order 
parameter profile are given in Subsec. iHBl followed by a de- 
tailed investigation of mean field results for the corresponding 
universal scaling functions in Subsec. llTCl In Subsec. iHDl we 
study the excess adsorption at corrugated surfaces. SectionHHI 
discusses critical Casimir forces between geometrically struc- 
tured walls mediated by the enclosed critical fluids. In partic- 
ular, Subsec. 1111 B I addresses the normal force between a fiat 
and a geometrically structured substrate, and in Subsec. UlI CI 
we discuss the lateral and the normal critical Casimir force be- 
tween two identically structured substrates. Section [TV] sum- 
marizes our findings. 



II. CRITICAL ADSORPTION ON GEOMETRICALLY 
STRUCTURED SUBSTRATES 

A. Model for geometrically structured substrates 

We model a laterally corrugated substrate confining a fluid 
by a periodic array of symmetric wedges and ridges as shown 
in Fig.Q] The direction along which the structure varies peri- 
odically is denoted as x, the direction indicating the distance 
from the substrate is called z, and the remaining (d — 2) di- 
rections in which the system is translationally invariant are 
denoted as y. Since for this system the observables do not de- 
pend on y, we are left with an effective two-dimensional prob- 
lem. The geometric structure of the confinement is character- 
ized (Fig. [TJi by the corrugation height h and the wedge angle 
y. The corrugation wavelength is given by l pro j = 2/ztan(y/2) 
and the actual substrate surface area is proportional to l S ide = 
2h/cos(y/2). The projected area A of the substrates in the 
(d — 1) directions {x,y} is macroscopically large (the nota- 
tion A — > °o will be suppressed in the following). 

One has to stress that all universal, system-independent 
quantities or functions discussed below reliably describe ac- 
tual physical quantities only for distances from the walls 
larger than typical molecular length scales of the system, 
i.e., the length scales of the constituent particles of the fluid 
and of the confining substrates. This means that all relevant 
length scales must be much larger than a typical molecular 
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FIG. 1: Side view on a periodic array of wedges and ridges modeling 
a geometrically structured substrate. The structure is periodic along 
the x direction and the shaded region indicates one unit cell. The 
direction globally perpendicular to the substrate is denoted as z. The 
structure is infinitely extended within the x-y plane, where y are the 
directions in which the structure is translationally invariant. 



size. For the correlation length | this requirement is fulfilled 
near criticality. For the geometric structure of the substrate 
this requirement implies that both h and also l pro j must be 
much larger than a characteristic, system-dependent molec- 
ular length scale. One should keep in mind that the require- 
ment that l pm j is larger than molecular sizes limits the range of 
validity of the universal properties for critical adsorption for 
small wedge angles y — > 0. 



B. General scaling properties for the order parameter profile 

In this subsection critical adsorption on a single geomet- 
rically structured substrate is studied. The order parameter 
<j> is a function of the spatial coordinates x and z, the re- 
duced temperature t, as well as of the corrugation parameters: 
= (p(x,z,h,y,t). Close to the bulk critical temperature (j) 
takes the scaling for m 14 ' 15 ' 47 ' 63 



<p(x,z,h,y,t) = a\t\ p P±(x±,z±,h±,y), 



(1) 



where x± —x/£,± is the direction of periodicity, z± = z/§± is 
the global normal distance from the substrate, and h± = h/B,± 
is the corrugation height. All lengths are scaled by the corre- 
lation length. The scaling functions P± are universal once the 
nonuniversal amplitudes a and §q are fixed as the amplitude 
of the bulk order parameter and of the true correlation length 
corresponding to an exponential decay of the pair correlation 
function. For distances from the substrate which are small 
compared to |±, or for T — ► T c , the scaling functions and the 
order parameter exhibit power law singularities 4 ^ 4 ^: 



<p{x,z,h,y,t = 0) =aC±(x/h,z/h,y) 



i) 



-p/v 



(2) 



where C± are dimensionless, universal amplitude func- 
tions depending only on distances scaled by the corruga- 
tion height h. The usage of the subscripts '±' for the 
scaling function in Eq. (O refers to the fact that they are 
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derived from the scaling functions P± in the limit t — > 
from above or below T c . Equation (f2| implies the relation 
C+{x/h,z/h, y) /C-(x/h,z/h, y) = (^ + /§T)~ fl/v between the 
scaling functions C± and the universal amplitude ratio of the 
correlation lengths above and below j^. 46 - 64 - 65 

1. Distant behavior 

For distances from the substrate much larger than the cor- 
relation length the scaling functions for T ^ T c decay as 

P+{x+,z+ -K*>,h + ,Y) « exp{-z+} (3) 

and 

P-(x-,z- ^ °°,h-,y) -1 « exp{-z_}. (4) 

Generally, for z/h — > °° it is not possible to detect the cor- 
rugation of the confining wall, and the system resembles a 
semi-infinite system with an effective planar wall. The lead- 
ing correction due to the corrugation is that, viewed from a 
distance, the wall does not appear as a planar wall located at 
z = but located at an effective height z = h e ff < h with 



z - h efj 



-1-/3 /v 



h eff=hheff(y), 

so that 

<j>(x,z,h,y,t = 0) =aC±(x/h,z/h,y) 
where 

C±(x/h, z/h, y) = C±(x/h,z/h,y) M -- fc e# (y) 



(5) 



(6) 



(7) 



The effective height is implicitly defined by the request that 
far from the wall 



C- 



0) y 



c±, 



(8) 



where c± are the universal surface amplitudes for the cor- 
responding semi-infinite planar wall system, defined via 15 

<j)°°/ 2 (z,t = 0) = ac± (z/^q,) , where the superscript 
°°/2 refers to the semi-infinite system. Off criticality the 
exponential decay of the scaling functions is proportional 
to exp{— (z— heff)/%±} so that even in the limit z — > °° the 
amplitude of the exponential decay °< exp(— z/%±) still car- 
ries the information about the corrugation via the prefac- 
tor exp(/z h e ff(y)/£,±). At T c the leading decay behavior is 
ac±(z/Eft)~P' v , which is independent of the corrugation, 
which for z — ► °° leaves its trace only through subdominant 
algebraic terms. 



C. Order parameter profiles within mean field theory 

Within a field theoretical renormalization group approach 
the general standard fixed point Hamiltonian for critical phe- 
nomena in confinements is given by 2 ' 3 ' 66 



where the dimensionless bulk Hamiltonian reads 

mm= ^A-|i(v^) 2 +^ 2 +^ 4 |. do) 

Equation ( TTOb holds for zero bulk fields as assumed in the fol- 
lowing. The integration in Eq. dTOb runs over the volume V 
which is accessible to the fluid, and r = {x,z,y} is the spatial 
position vector. In Eq. (O the surface Hamiltonian Jf? s and the 
edge Hamiltonian J%? e incorporate the surface (edge) enhance- 
ment of the coupling energy at the boundary and the surface 
(edge) fields which act on particles close to the surface (edge) 
only. At the critical adsorption fixed point 5 the surface and 
edge contributions turn into boundary conditions correspond- 
ing to infinite surface and edge fields describing the strong ad- 
sorption limit so that 1 surface = +°° at all surfaces and edges 
of the system. Since our study deals with identical chemical 
boundary conditions at all surface and edges, it is sufficient 
to consider only the '+' case corresponding to (j) = +°°. In 
Eq. ( [Tol l the coefficient T is proportional to the reduced tem- 
perature t, and the coupling constant u is positive. 

Quantitative results can be obtained within the framework 
of mean field theory (MFT), which corresponds to the zeroth 
order contribution in a systematic £ = A — d expansion within 
the field theoretical renormalization group approach. Within 
the MFT one neglects fluctuations of the order parameter, and 
only the configuration m — \J (h/3!)(0) with the largest sta- 
tistical weight exp(— [0]) is considered. This configuration 
fulfills the Euler-Lagrange equation 



V m = Tin 



(11) 



which is obtained via functional minimization of Eq. ( TTOb . 
The universal quantities and scaling functions calculated via 
MFT are exact for dimensions larger than the upper criti- 
cal dimension d uc — 4. The mean field values of the crit- 
ical exponents j3 and V are j3 = V = \, and for the coef- 
ficient T in Eq. (TTOt one has within MFT T = t (^ + ) -2 for 
t > and T = 2?(§o~) -2 for t < with the universal ratio 
(^+/^) 2 = Within MFT and d = 4 the coupling con- 

stant u in Eq. ( TTOb c an be re lated to the nonuniversal constant a 
in Eq. ([T]) as a = y/ (3!/k)/^q". For T = T c the dimensionless 
value m x h depends on x/h, z/h, and y only and is indepen- 
dent of the nonuniversal quantities a and Ej" (see Eq. (0). 

The order parameter profiles presented in the following 
have been calculated numerically. In order to obtain a bound- 
ary condition for the numerical calculation we use the short 
distance expansion of the solution of Eq. (fTTb for the semi- 
infinite system with the boundary condition m = +°° at the 
wall 



V2_ 
R 3\/2 



R + 0(R 3 ), 



(12) 



j^[<l>] = je b [<j>}+j^[<j>]+jf e [(j>}, 



(9) 



where R is the minimal distance to the wall. Eq. (fT2b holds 
only for distances R -C ^,h,l pru j. We approximate edges by 
small inscribed circles and use Eq. (fT2l) for the numerical 
boundary condition close to edges with the distance R mea- 
suring the minimal distance to the edge. We find the universal 
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FIG. 2: Order parameter scaling function P+ (x+ , z+ , h+ , y) for 
T > T c , where x + = x/%+, z+ = an ^ ^+ = V?+> m terms of 

contour plots for geometries with constant scaled corrugation height 
/j+ = 3 but varying wedge angle y = 19° in (a), 90° in (b), and 128° in 
(c). For small values of y and for increasing normal distances z+ the 
profile rapidly adopts an effective planar wall behavior correspond- 
ing to straight contour lines. For large values of y the wall structure 
reaches deeply into the bulk. 



scaling functions for the order parameter profiles upon func- 
tional minimization by using a conjugate gradient algorithm 
within a finite element method. The grids and the numeri- 
cal boundaries together with the minimization conditions have 
been chosen such that for the order parameter profile the rela- 
tive error is less than 0.1%. 



1. Comparison of order parameter profiles 

In terms of the universal scaling function P + (see Eq. ([1)) 
examples of order parameter profiles for T > T c are shown in 
Fig. |2 (In the case T < T c , the scaling function P — 1 looks 
very similar provided /z_ = h + and y are the same.) For con- 
stant values of h± the influence of the corrugation on the order 
parameter profile along the z direction varies strongly with the 
wedge angle y. For small y, upon increasing z+ the profiles 
rapidly adopt an effective planar wall behavior correspond- 
ing to straight contour lines of the order parameter. For large 
wedge angles y the corrugation induces a spatial variation of 
the order parameter also far from the wall. That is, the depen- 
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FIG. 3: Comparison of the contour lines of the scaling function P + 
for T > T c (a) and of the order parameter profile m for T = T C (b), 
respectively, for the same wedge angle y = 90° but for various cor- 
rugation heights h. The values of P+ in (a) and m in (b), respectively, 
have been chosen such that all contour lines coincide at x/h = 1. 
(The values of P + in (a) are 1.0 for h/%+ = 1, 0.093 for h/t,+ = 3, 
and 0.0005 for h/£+ = 8; the values of m in (b) are 1.61 x 10" 2 /<^ + 

forVI + = 100 ' 53 x 1(r3 /£ + for h l%a = 300 ' and2 - 0x 10-3 /<? + 
for h/^Q = 800.) For fixed y, in (a) the order parameter scaling func- 
tion for T > T c depends on the values of z/|+, */£,+ , and h/^+ (or 
equivalently on z/h, x/h, and in contrast to the profile for 

T = T c in (b), which depends only on the distances z and x scaled by 
the corrugation height h, so that for various values of h the contour 
lines coincide, which amounts to a welcome check of the numeri- 
cal data. (The slight deviations of the contour lines in (b) are due to 
numerical errors in the determination of the contour lines.) 



dence of the order parameter profile on y is crucial and cannot 
be split off in a simple way. 

A comparison of contour lines for the order parameter pro- 
file at a constant wedge angle y but for varying corrugation 
height is shown in Fig. [3]for T > T c and T = T c . For T = T c 
the order parameter depends only on the scaled lengths x/h 
and z/h [Fig. [3jb) and Eq. ©J in contrast to the off-critical 
case, for which the order parameter depends in addition on 
h± = h/E,± [Fig.|3a) and Eq. (Q}]. For increasing values of 
h + =hj%+ the undulation of the contour lines does not only 
increase in units of z/^+ but also relatively in units of z/h as 
can be inferred from Fig. IHa). In Fig. |4]a comparison of the 
order parameter profile for the corrugated wall geometry with 
the analytic expression 46 for the order parameter profile close 
to a single wedge/ridge is shown. As expected, close to the 
edges the analytic solution coincides with the results for the 
geometry of a periodic array of wedges and ridges. 
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FIG. 4: Contour plot of the order parameter profile m at criticality for 
y = 90°. The gray solid lines correspond to values of m x h = 1.6, 
2.4, 3.2, 4.8, 8.0, 16, and 48 (top down). For T = T c the value of 
in x h depends only on x/h, z/h, and y. The corresponding analytic 
solutions 46 for single wedges and ridges are shown as dashed lines. 
Close to the edges the order parameter profile for a single wedge or 
ridge coincides with the one close to a periodic array of wedges and 
ridges. 
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FIG. 5: Sketch of the undulation As of the contour lines of the or- 
der parameter profile measured in units of the corrugation height h. 
As depends on the value s = z/h of the mean position of a contour 
line along the z direction. The undulation of the contour lines is a 
measure of how deeply the surface structure reaches into the bulk, 
independent of the absolute value of the order parameter itself. For 
s — > oo the undulation As vanishes. 



2. Proliferation of the surface structure into the bulk 

In order to discuss the propagation of the substrate structure 
into the bulk in more detail, we study the behavior of the con- 
tour lines of the order parameter. The deviation of a contour 
line from a straight line gives a measure of how the structure 
influences the adsorption behavior at some distance z from the 
surface, independent of the actual value of the order parameter 
at this position. Using this measure, the cases T < T c , T = T c , 
and T > T c can be compared in a straightforward way. To this 
end we introduce, in units of h, the undulation As(s) being the 
width of a given contour line along the z-direction as a func- 
tion of its corresponding mean position s = z/h (see Fig. |5J. 

As a distant wall approximation we introduce a superposi- 
tion ansatz for the order parameter, 

m(x,z,h,Y,t) ~m(°\z,t) +m( 1 \x,z,h,y,t), (13) 



where is the effective planar wall solution, and is the 
corrugation contribution to the order parameter profile, which 
is assumed to factorize into x- and z-dependent part as 



2 (1) (x,z,/!,y,f) = m (jc 



(x,h,Y,t)xm^(z,h,Y,t). 



(14) 



Inserting the ansatz in Eq. (fT4l into the Euler-Lagrange equa- 
tion [Eq. ([TTl il and using the fact that far from the wall, 
where s> 1 and As(s) <§; 1, the corrugation contribution is 
much smaller than the planar wall contribution, one finds for 
the major corrugation contribution a sinusoidal dependence 
mW oc cos(2tzx/ Iproj) along the x direction, and that m$ is 
decaying exponentially along the z direction for all temper- 
atures, including T = T C . This is discussed in detail in Ap- 
pendix [A] For T = T c this is in contrast to the planar wall 
contribution, which at T c decays algebraically. For s ^> 1, i.e., 
within the distant wall approximation, the undulation of the 
contour lines behaves as (see Appendix lAl or Ref.l67l) 



As(s) 



(s-Mr)) exp{-CT (y)i}, 
exp{-a±(V£fc,y)*}) 



T = T c , 
T § T e , 



(15) 



where Ob and <j± characterize the exponential decays: 

(7 (y) - 7rcot(y/2) (16) 

and 



o±(A±,y) =h± 



1 + 1 — cot(y/2)) -1 



(17) 



These expressions are obtained by assuming that the angle y 
is not too small and, in the case t ^ 0, that h and l pm j are 
comparable with £±. 

The decay constant Ob(y) for t = corresponds to an ex- 
ponential decay oc exp{—2nz/l pr oj} of the amplitude of the 
order parameter undulations. This decay is analogous to sub- 
strate potentials near crystal surfaces with a periodic atomic 
structure, where that part of the potential which is due to the 
atomic corrugation decays exponentially, too, although the 
strength of the substrate potential decays algebraically (see, 
e.g.,Refs.|68||69l). 

In addition, the undulation of the contour lines can be de- 
termined from numerically calculated order parameter pro- 
files. The numerical data confirm that the undulations decay 
exponentially for all temperatures. Moreover, Fig. [6] shows 
that the numerically determined decay constants compare well 
with the results obtained within the distant wall approxima- 
tion (Eqs. ( fToT i and (fTTIi); there are significant deviations only 
outside the range of validity of the approximations used to de- 
rive Eqs. ( fT6] l and ( fT7] i. Interestingly, the undulations of the 
contour lines for T = T c decay more rapidly than for T ^ T c , 
in contrast to the order parameter profile itself. Furthermore, 
as mentioned above, Eq. dTvT > implies that for large values 
of h / £± the effect of structuring the substrate propagates not 
only in absolute units but also in units of h deeper into the 
bulk than for small values of hj%±. 
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FIG. 6: Decay constants do (T = T c ) and a± (T ^ T c ) correspond- 
ing to the exponential decay of the undulation As(s) (see Fig. O of 
the order parameter profiles for critical adsorption on a single pat- 
terned substrate (see Fig. [TJ. The analytic distant wall approxima- 
tions (lines; Eqs. U6\ and U7» compare well with the full numerical 
data for the order parameter profiles (symbols). There are devia- 
tions outside the range of validity of the underlying distant wall ap- 
proximation. Within the numerical accuracy the analytic prediction 
a + (h + ) = d-{h ) (Eq. (XT}) is confirmed. 



For certain substrates it is possible to keep 1^, i.e., the 
actual substrate surface, constant, but to vary the wedge angle 
y, like an accordion. Experimentally, such tunable, periodic 
buckled surfaces can be created using wrinkling polymer films 
under stress i 70 ' 71 For the periodic array of wedges and ridges, 
the unsealed undulation of the contour lines Az(z) = hAs(z/h) 
is maximal for y max fulfilling 



dAz(z) 
dy 



= 0. 



(18) 



7= 7m. 



Within the distant wall approximation given above and for 
T = T c , upon neglecting all prefactors of the exponential de- 
cay of As Eq. ( fT8b becomes 



sin 3 (y m , x /2) _ 2nz 
cos 2 (y„ MX /2) I 



(19) 



side 



This means that the angle y max for maximal undulation de- 
pends on z/hide, and for z/hide — ► 00 it attains n, i.e., if one 
wants to achieve significant undulations even far away from 
the substrate its wedges have to be wide open. 



3. Distant behavior 

As shown in Fig. [7] for T = T c the order parameter pro- 
file takes the effective planar wall form for z/h 3> 1. This 
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FIG. 7: Universal scaling function C±(x/h,z/h,y) for the order pa- 
rameter at T = T c [see Eq. Q] in units of the planar wall universal 
surface amplitude c± [see Eq. (8)]. The lines represent order param- 
eter profiles corresponding to the lateral position of the wedge center 
(x = 0) and the edge crest (x = l pro j/2) for various values of y [see 
Fig-UJ. Far from the corrugated substrate (z/h ^> 1) the order pa- 
rameter resembles the one close to a planar wall located the position 
z = h e ff. The effective substrate height h^g is determined such that 
C±/c± — > 1 for large z/h [see Eqs. (|5}h[8}]. The inset shows the di- 
mensionless effective height h e ff(y) = h e ff/h (symbols). The dashed 
straight line 1 — y/(2n) approximates the data well. 



means that far from the corrugated wall the order parameter 
resembles the one for a planar wall located at z = heff. Based 

on Eqs. (O and © the effective height h e ff(y) introduced in 
Eq. © is determined from the order parameter profiles. From 
the numerically determined data shown in the inset of Fig. [7] 
we find that 



(20) 



For y — * the effective height is h e jf = h, which means that the 
'closed' wedges form a planar wall located at z = h. For y — * K 
the effective height attains the value of h e g = h/2, and the dis- 
tant profile resembles the one created by a planar wall located 
at the mean corrugation height. These limits are consistent 
with results obtained for the quantum electrodynamic Casimir 
force between corrugated surfaces. 51 In Fig. [8] a cut through 
the order parameter profile for T ^ T c is shown. The order pa- 
rameter profile exhibits an effective planar wall behavior with 
an exponential decay towards its bulk value forz ^ h. It turns 
out that the effective height for a given wall structure is inde- 
pendent of temperature, i.e., the same for T = T c and T ^ T c , 
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FIG. 8: Variation of the order parameter (Eq. dl) along the z direc- 
tion near a periodic array of wedges and ridges (see Fig.[TJ for T >T C 
(a) and T < T c (b). The corrugation is chosen as h± = 3 and 7 = 90°. 
The dashed lines correspond to order parameter profiles at the lat- 
eral position x = 0, i.e., at the center of the wedge, the solid lines 
correspond to the lateral position x = l pro i/2, i.e., at the crest of the 
edge (see Fig. [TJ. For z/h^> 1 these profiles join, which means that 
far from the surface the order parameter exhibits an effective planar 
wall behavior. The order parameter decays exponentially towards its 
bulk value P + (x + , z+ — » °°,h + ,y) = 0, or P-(x-,z- — > °°,/?_ , y) = 1, 
respectively. The gray lines correspond to the analytic solutions of 
the order parameter close to a planar wall 2 '- located at the position 
z = h e g = hh e ff(Y) (dotted vertical lines) where they diverge. The 
effective substrate height h e jfis determined according to Eq. d20b . It 
turns out that the value of hgte for a certain wall geometry is the same 
for T = T C and T 7^ T c . 




■K °< exp{— z.-} 

x- =0 

X_ = l pm j/ {2$-) 

coth{{z--hJi eff {y))/2}-\ 



z smaller than or comparable with a typical molecular length 
ao. However, for larger distances z > «o one hasi£ 



a 

^ /2 (o= Jdz{r /2 (z, 



■Hz-- 



- 00 t 



'} 





^ p - v j dz ± {pZ'\z±)-PZ' 2 H}- (22) 



On the molecular scale the order parameter is bounded for 
z — ► 0, whereas the divergence of the universal scaling func- 
tions reflects the continuum description. Thus, the first inte- 
gral is a finite nonuniversal quantity, which is subdominant to 
the second, diverging integral. The universal amplitude of the 
singular part of the excess adsorption of a planar wall is 



r+ /2 = lim 



0/2 



(O/(0(z = «v) 4±) 



'dz±{p± /2 (z±)-P± /2 H], (23) 



where the lower limit of the integral has attained zero for 

For a periodic array of wedges and ridges it is sufficient 
to consider a single unit cell as sketched in Fig. Q] For one 
unit cell the excess adsorption T± per [d — 2) dimensional 
projected area is given by 



r±(h,Y ) t) = a(tf) 2 \t\i } - 2v r±(h ±) y), 



(24) 



with its scaling function 



dx± I dz± 



r±(A±,y) = 

{p± (je± ,z± , A± , 7) - p± /2 (25) 
For a planar substrate, i.e., 7 = 7T, Eq. d25l l reduces to 



f ± (/, ± = 0,7=7r) = ^fl/ 2 . 



(26) 



D. Excess adsorption 

Certain experimental approaches such as adsorption mea- 
surements provide access only to integrated quantities rather 
than spatially resolved local quantities. For a binary liquid 
mixture of A and B particles close to its critical demixing 
point at T c the overall enrichment of, e.g., the A particles com- 
pared to the B particles at a planar wall with perpendicular 
direction z is given by the excess adsorption per unit wall area 



«/2 



(0 



£dz{r /2 (z,t)-<j>(z^^,t)}. 



(21) 



where the subscripts ± refer to T ^ T c , and the superscript 
oo/2 refers to quantities in a semi-infinite system. The scal- 
ing form of the order parameter does not hold for distances 



In order to compare the critical adsorption on geometrically 
structured substrates with the critical adsorption on planar 
substrates, we define a reduced relative excess adsorption, 



r±(ft ± ,y)-^r; 
1 fW2 



/2 



(27) 



T± relates the excess adsorption for a unit cell of a corrugated 
substrate to the excess adsorption for a planar substrate of lat- 
eral extension / in x-direction in units of the latter one. Ob- 
viously, there are two interesting choices for I: the projected 
width of a periodic wedge l pro j, and the actual surface length 
I side °f a un it wedge. 

For these limiting cases one can obtain explicit expressions 
for T±. For h± = h/%± — > 00 the corrugation resembles a 



8 



collection of large tilted planar walls of length / v( y e with a 
relatively thin layer within which the order parameter devi- 
ates from its bulk value. For large l s [ C \ e the edge contributions 
within a unit cell become unimportant. Likewise, in the limit 
7 — > % the wavelength of the corrugation is very large and the 
edge contributions can be neglected, too. Thus, in those lim- 
its the adsorption in one unit cell of the array of wedges and 
ridges is the same as the adsorption at a planar wall of lateral 
size l S ide, and one has 



h = - 



or y-^Jt 



hide lproj 
lproj 



1 



sin 



-1. 



(28) 



On the other hand for h± — > the corrugation becomes van- 
ishingly small on the scale of the correlation length, and the 
total adsorption in a unit cell is the same as the one of a planar 



wall of lateral size l pro j. Therefore, one finds 



l m\ h± ^°> lproi ~ lside 

5 > hide 



(I)-l- (29) 



For 7^0 and arbitrary values of h± we have not found a 
simple asymptotic formula. 

With Eq. ( l27i i we define the ratio R of the excess adsorption 
above and below T c as 



R(h+,y) = 



r+(h + ,r) 



r 



(h- = ft H 

/2 



• r) 

i+t 4 



(30) 



Mi) 



^ f: 



'12 



l+T-fA- 



which is independent of the choice of I. r~ /T_ is a uni- 
versal constant for critical adsorption at a planar wall and is 
experimentally accessible.— 



1. Universal behavior of the excess adsorption within mean field 
theory 

In the remainder of this subsection we present the results 
for the excess adsorption obtained from order parameter pro- 
files calculated within MFT Due to the divergence of the or- 
der parameter at the surfaces in the limit of strong adsorption, 
which in this scaling limit leads to a divergence of the excess 
adsorption, one has to calculate T± carefully. T± is obtained 
by introducing a cut-off for small distances from the substrate 
and subsequent extrapolation for the cut-off sent to zero. 

Figure|9]displays the reduced relative excess adsorption for 
t > in comparison to a planar wall of lateral size I — hide', 
in Fig. [TO] the comparison with a planar wall of lateral size 
I = lp ro j is shown. The general asymptotic behaviors derived 
above (Eqs. (|28l l and d29l i) are attained in all cases. Concern- 
ing the interesting behavior for small wedge angle 7, the to- 
tal excess adsorption is less than the one for planar walls of 
lateral size l s ^ e . This means that for all corrugation parame- 
ters h± and 7 the loss in adsorption near an edge exceeds the 
gain in adsorption inside a wedge 47 ; this effect becomes more 
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FIG. 9: Reduced relative excess adsorption T+ within MFT at a sin- 
gle geometrically structured substrate [Fig.Q") for T > T c [Eq. J27H. 
comparing a structured geometry with a planar wall of lateral size 
hide- Note that l s ^ e = 2h/ cos (7/2). The lines connecting data points 
are smooth fits, while the thick dashed-double-dotted line shows the 
asymptotic behavior for h + — * [Eq. {29)]. The negative range of 
values of Y + means that the amount of adsorbed matter on a geomet- 
rically structured substrate is less than the one on a planar wall with 
the same actual substrate surface size. All curves meet at the point 
{y = 0,Y+(h+,Y,l s ide/^+) = ~ 1}i because for 'closed' wedges the 
excess adsorption in a unit cell of the periodic array of wedges and 
ridges vanishes compared with the reference value of the excess ad- 
sorption on a planar substrate with the same actual surface size. 
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FIG. 10: Reduced relative excess adsorption T + within MFT at 
a corrugated substrate [Fig. [JJ for T >T C [Eq. d27H in compari- 
son with the one on a flat surface of lateral size l pm j- Note that 
lp ro j = 2/i tan(y/2). The lines connecting data points are smoothly 
fitted, whereas the thick dashed-double-dotted line shows the asymp- 
totic behavior for 7 — > 71 or h + — > 00 [Eq. J28H. The excess adsorp- 
tion on a geometrically structured substrate is larger than the one on 
a planar substrate of the same projected area. For 7 — > the ratio 
of these two quantities even diverges. This is in contrast to the case 
/ = l side in Fig.[9] 
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pronounced upon approaching T c , However, for y — > the 
amount of adsorbed matter is larger than at a wall of lateral 
size l pr0 j. Indeed, Y±(h/%±, Yjpmj/^±) even seems to diverge 
for y^O, i.e., the adsorption on a geometrically structured 
substrate with a short corrugation wavelength is much larger 
than the one for a planar substrate with the same projected 
area. This effect becomes more pronounced away from T c . 
Thus, if one aims at enhancing the effect of critical adsorp- 
tion, one should structure the adsorbing substrate. At a first 
glance this appears to contradict the conclusion above that the 
order parameter profile for z > h and for y ^ closely resem- 
bles that of a planar wall [Fig.[2J a )]- However, for T ^T c and 
within MFT the major contribution to the adsorption stems 
from the vicinity of the substrate. This part diverges due to the 
increase °= R 1 of the order parameter profile near the surface, 
where R is the closest distance to the surface. That is, within 
MFT the major contribution to the adsorption is governed by 
the wedge walls. Since the wedge walls are accessible for the 
fluid even for y ^ and because the size of the wedge walls 
is larger than l pro j, the relative excess adsorption diverges for 
y — > (which implies l pm j — ► 0). However, in order that the ac- 
tual physical quantities are described properly by the universal 
scaling functions, l pro j must be large compared with molecu- 
lar length scales, which imposes a limit on the smallest wedge 
angle below which the universal features do no longer capture 
the structure near the wedge center. 

Beyond MFT, the behavior of Y±(h±,Y,lproj/%±) is ex- 
pected to be qualitatively different. In this case the critical 
exponent ratio for the order parameter decay at T c and at a 
planar substrate (Eq. ©) is j3/v ~ 0.5 (instead of j8/v = 1 
within MFT), and the divergence of the excess adsorption 
stems from the slow decay far from the surface and occurs 
only for T — ► T c . Beyond MFT, there is no diverging contri- 
bution from the vicinity of the substrate. Since for small y the 
order parameter profile outside the wedges of the periodic ar- 
ray of wedges and ridges resembles the one of a planar wall 
(see Fig.|2a)), one expects the excess adsorption in a unit cell 
of the periodic array of wedges and ridges to approach for 
y — > the value of the excess adsorption at a planar wall of 
lateral size l pm j. That is, Y±(h±,y — > 0,l pro j/^±) — > beyond 
MFT (in contrast to the divergence within MFT [Fig. [TOl ). 
However, for intermediate wedge angle ranges < y < n the 
excess adsorption at a corrugated substrate is larger than at a 
flat one, i.e., generally Y±(h±, Y,lproj/£,±) > 0. On the other 
hand the behavior of Y±(h±,y, hide/ %±), which measures the 
difference of the excess adsorption at a corrugated substrate 
and a planar wall of lateral extension hide, should be quali- 
tatively similar to its behavior within MFT [Fig. |9) because 
also within MFT a major contribution to this difference is 
due to behavior distant from the walls. Beyond MFT, due 
to the slow decay of the order parameter far from the wall for 
T — > T c , the property of the excess adsorption at a corrugated 
wall to be smaller than the corresponding one at a planar wall 
of lateral extension hide is even more pronounced than within 
MFT. That is, also beyond MFT T± (A± , 7, hide /£±) < 0, and 
Y±(/j±,7^0,W§±)^-1. 

For t < the reduced relative excess adsorption looks sim- 
ilar to the case t > 0. However, the curves for h- — h + do 
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FIG. 11: Universal ratio R of the excess adsorptions at a single 
topologically structured substrate [see Fig. Q~) above and below T c 
[Eq. j30H, divided by the corresponding planar wall ratio and /E, + . 
The scatter of the data is due to the difficulty to interpolate between 
the data points. Besides these uncertainties the curves for different 
h + = A_ and / = l s id e >lpmj do not differ much. Approximately, for all 
parameters the data vary as a function of y as ~ 1 — 0.07[(y— n) /it] 4 
(shaded region). The ratio of the excess adsorptions above and below 
T c is smaller at rough substrates than on planar ones, especially for 
small values of 7. 



not fall on top of each other. One can determine the ratio R 
of the excess adsorption above and below T c as introduced in 
Eq. (f30b by studying these differences. The dependence of R 
on various parameters is shown in Fig. QT| According to its 
definition in Eq. ( f3Qb . R is independent of the choice for /. In 
addition, it seems to be independent of the choice for h+, too. 

5- 



R attains 



=±77 for 7- 
rj 



it, but decreases for 7^0. Thus, the 



ratio of the excess adsorption on corrugated substrates above 
and below T c is different from the corresponding planar wall 
ratio. 



III. CRITICAL CASIMIR FORCES BETWEEN 
GEOMETRICALLY STRUCTURED SUBSTRATES 

In this section we discuss critical Casimir forces mediated 
by critical fluids confined between two opposing geometri- 
cally structured walls. Two types of systems are considered: 
First, we study the normal critical Casimir force between a 
geometrically structured wall and a planar wall. Second, we 
discuss the lateral and the normal critical Casimir force be- 
tween two identically structured walls. Throughout this sec- 
tion the chemical boundary conditions are the same for the 
two walls, i.e., the signs of the surface fields at the two sub- 
strates are equal. As before, the geometric structure of the 
substrates is taken to be a periodic array of wedges and ridges 
along the same lateral direction x. The corrugations are on top 
of walls separated by a distance L in z-direction (see Figs. [12] 
and, c.f., [I6l l, As before, the projected (d — 1) dimensional 
area A of the two substrates within the {x, y} plane is macro- 
scopically large. If both substrates are structured, the wedges 
and ridges of the two substrates can be shifted by D relative to 
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each other. Thus, both cases, a structured substrate opposing 
a planar one, or two identically structured substrates opposing 
each other, are described by the variables L, h, and y (and D 
for the latter case), in addition to the reduced temperature t as 
a thermodynamic variable. 



A. Stress tensor in periodically structured confinements 

The free energy of such a confined fluid system decomposes 
into four distinct contributions^: 



F = VF b + SF s +EF e + 8F, 



(31) 



where V is the volume accessible to the fluid, S denotes the 
actual substrate surface area, E is the sum of the edge lengths, 
and F/, s e are the bulk, surface, and edge free energy densities. 
8F corresponds to that part of the free energy stemming from 
the finite size effect and the effect of corrugation. In the fol- 
lowing we focus on the singular part of F near the bulk critical 
point (without introducing a separate notation for this part of 

The normal critical Casimir force acting on the confining 
walls is defined as 



f± = 



1 



k B TA dL 



SF, 



(32) 



which does not include the bulk contribution to the force 
— kgTA 7E 7 wr, ich is constant with respect to L, and the lat- 
eral critical Casimir force is 



f\\ 



1 



k B TA dD 



SF. 



(33) 



By using the stress tensor, the forces can be calculated di- 
rectly from the order parameter profiles. This has the advan- 
tage that one does not face the numerical difficulties of cal- 
culating differences of free energies which diverge due to the 
divergence of the order parameter profiles in the scaling limit 
near the surface. For the fixed point Hamiltonian given by 
Eqs. (0 and (TlOb the stress tensor components ar o 10 i 73 ' 74 - 75 



T U (t) 



dip d<j) _ 

"5 "5 °kJ 

dr^ ari 



-4/(r), 
(34) 

where k,l = x,z, ••■ indicate the components of the position 
vector r = {x,z, ■ ■ •}, and 4/ is the improvement term, which 
reads 



4z(r) 



d-2 
4(d- 1) 



2a2 



d 2 <j> 

dnd 



r k dr, 



(35) 



For the systems under consideration, T zz and T xz correspond to 
the normal and lateral critical Casimir forces, respectively, as 
introduced above if bulk contributions are subtracted (within 
MFT these occur only for T < T c ). 

In the following we only consider confinements with over- 
all periodicity. This is the case for the corrugation model used 
here with an overall periodicity wavelength l pro j if the ratio 
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FIG. 12: Confinement composed of a structured wall and a/lat wall 
(denoted as (sf) geometry). The wedge vertices of the geometrically 
structured substrate are separated by the distance L from the planar 
wall. 



of the two single corrugation wavelengths is a rational num- 
ber or if one of the two substrates is fiat. Then the force in 
fe-direction per unit area is 



f 1 

/* — ; — 

lproj 



(36) 



where is evaluated at a fixed position z = Zo in between the 
two substrates. Due to the periodicity of the system, the con- 
tribution to the force stemming from the improvement term 
[Eq. ((35])] vanishes 67 : 



1 fh'ivj 
- — / dxlfo = 0. 



(37) 



Normal critical Casimir force between a geometrically 
structured and a flat substrate 



As a first step for studying the effect of corrugation on criti- 
cal Casimir forces, this subsection deals with a structured wall 
opposing a/lat wall as shown in Fig. [HI denoted as (sf) ge- 
ometry. Within the (sf) geometry there is no lateral force. 
The normal critical Casimir force takes the following scaling 
form: 



f ( l f) (L,h,y,t) 



(d-l)L- d fl sf 2(y ± ,h ± , Y ) 



(38) 



(d-l)IT d A^ 



:(*/) 



where y± —t(L/£,^y^ v , and/i± =h/t;±; + is the universal 

""(sf) 

scaling function of the normal force, A ± is the universal gen- 
eralized Casimir amplitude, and A+,+ is the universal Casimir 
amplitude for two planar plates with parallel surface fields.^i 
For MFT d = 4, and y+ = xL 2 for T > T c , and y_ = 2 xL 2 for 
T < T c ; the Casimir amplitude is A+.+ ~ -15.7561 (3!/m). 10 



For T = T c the dimensionless quantity f ± 
L/h and 7, only [Eq. (|38)]. 



sf) 



x h depends on 
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1, Distant behavior 

The large distance behavior of the normal force is governed 
by the effective planar wall behavior of the corrugated sub- 
strate. That is, the scaling functions of the normal critical 
Casimir force resemble the one for two planar walls at a dis- 
tance L e ff = L — heff (see Eq. d38b ): 



and 



fr±>(y±,h±,Y) 



(»/) ( L J\ i>>/! ' / ' 



ti) <39) 



1 - 



i -i 



-d 



(40) 



where y± ff = t(L eff /^) l l v and /+'+ (y^) is the scaling func- 
tion for the critical Casimir force between two planar walls 
at distance L ej y . 6 ' 7 ' 8 ' 9 ' 10 i 24 (In principle h e g = h h e ff{y) for 
the confined system could differ from the one introduced in 
Eq. (0 for the corresponding semi-infinite system; however, 
it turns out that they are the same.) 

In order to provide quantitative expressions for the univer- 
sal scaling functions we calculate the stress tensor compo- 
nents within MFT (see Subsec. lHCl l corresponding to the ze- 
roth order of the e-expansion. These numerical results for the 
scaling functions of the normal force are shown in Fig. [13] 
where the effective height is the same as the one for a single 
corrugated substrate [Eq. d20b l. The data agree very well with 
the effective planar wall limit for large values of L/£,±. We 
find that for constant h± the effective planar wall limit of the 
normal force is reached already for smaller values of L/%± if y 
is smaller (see Fig.[T3rb)). On the other hand, we find that for 
constant y the effective planar wall limit is reached faster as 
a function of L/£± if the value of h± is smaller. Quantitative 

results for the generalized Casimir amplitude A\j^ [Eq- d38l ll 
are presented in Fig. [14] The numerical data agree with the 
effective height limit (Eq. d40b together with the semi-infinite 
result for h e ff(y) in Eq. (|20| >) for large values of L/h. The ex- 
pression for the effective height [Eq. (f20b l is consistent with 
results for the quantum electrodynamic Casimir force between 
a corrugated and a planar wall. 51 



2. Nearby behavior 

As indicated in Fi g s . [T3[ and [T4l there are pronounced devi- 
ations from the effective planar wall limit if the wall distance 
is small, i.e., for L t i p = L h — ► 0. If the characteristic sizes 
of the corrugation are much larger than the distance L t i p , a 
Derjaguin-like approximation may be appropriate. 76 Within 
this Derjaguin-like approximation we replace the corrugation 
of wedges and ridges by a set of staircases with infinitely small 
horizontal terraces and vertical steps and sum the single con- 
tributions of the terraces in order to obtain the total normal 
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FIG. 13: (a) Universal scaling function of the normal critical Casimir 
force (in units of the absolute value of the universal critical Casimir 
amplitude A+.+ for planar walls) for the (sf) geometry [Fig. 1121 
within MFT, for constant h± = 1.5 and for three wedge angles y, as a 
function of y± =f(L/^) 1//y . The gray line corresponds to the MFT 
result for the normal critical Casimir force between two planar walls 
at distance L. The absolute value of the force between a geometri- 
cally structured substrate and a planar wall at distance L [Fig. 1121 is 
stronger compared to the one between two planar walls at distance L, 
because the corrugation on top of the lower substrate in Fig. [inef- 
fectively decreases the distance between the two opposing substrates. 

The geometric constraint L> h [Fig. 1121 implies \y±\ > }ij_ v and 
thus leads, for h± fixed, to a divergence of the scaling function of 
the normal critical Casimir force for \y±\ \ \yi,±\ = \t\(h/^y^ v , 
where the planar wall approaches the tips of the corrugation; for 
fixed h± smaller absolute values of y± cannot be reached for a given 
(sf) geometry. From Eq. ([43} one infers that, for fixed h± ^ 0, 

f±,± (y± "> Vh± , *±, r) « I [y±l - b± I \~ v{ - d - 1 \ On the other hand, 
for / — » both h± and y± vanish so that in this limit the scaling vari- 
able y± can reach zero, leading to the finite result for /£ ± given 
by the second line of Eq. ( 1381 ). (b) Rescaled (i.e., multiplied by 
(L e ff/L) ; see Eqs. ( 1391 ) and ( 140b ) universal scaling function of the 
normal critical Casimir force within MFT (d = 4) for the same con- 
finements as in (a) as a function of ±L e ff/t,±, where ± indicates 
T T c . Note that the curve for y = 37° diverges to — °°, too; however, 
due to numerical limitations, here only the indicated limited range is 
covered. The gray line corresponds to the normal critical Casimir 
force between two planar walls at distance L e g (see Fig. [12] and the 
main text); note that within MFT y^g = TL^, and y~g = 2tL^. For 
L/h ^> 1 the normal critical Casimir force between a corrugated and 
a flat substrate reaches this effective planar wall limit [Eqs. J39b and 
( I40H with L e ff = L — h e ff and h e g = hh e ff(y) given by the semi-infinite 
expression for h e g(y) [inset of Fig.[7j. 
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FIG. 14: Universal generalized Casimir amplitude A ± (L/h,y) 
[Eq. J38H of the normal critical Casimir force at criticality for the 
(sf) geometry (see Fig. 112b as a function of L/h. The results cor- 
respond to MFT (d = 4). In this representation the critical Casimir 
force between two planar walls at distance L corresponds to unity 
(horizontal line). For L/h^ 1 the normal critical Casimir force 
at T — T c approaches the corresponding effective planar wall limit 
[Eq. J40H with L e g= L—hh e ff(y) where h e ff(y) is given by the semi- 
infinite expression [inset in Fig. [7). The effective planar wall limits 
would diverge at L/h = h e ff(y) < 1 as (L/h — h e ff(y))~ d . However, 
this approximation becomes unreliable for small distances. The ac- 
tual scaling function diverges at L/h = 1 as (L/h — l) _ ( d_1 ) 
(see Eq. d44l>). For wide wedge angles y the deviations from the ef- 
fective planar wall limit are most pronounced. 



force: 



JDer) 



proj 



dx ■ 



(L(x)) 



(41) 



where L(x) = L t j p + h(l — 2\x\/l pm j) is the local wall distance 
and ft±(y±(x)) with y±{x) = {L{x)/t, ± ) l l v is the normal 
critical Casimir force scaling function between two planar 
walls at distance L{x). The summation of single contributions 
is clearly a strong approximation to the actual non-local char- 
acter of critical Casimir forces and it is reasonable only in 
the nearby range where local force contributions govern the 
behavior. Performing an integration by parts in Eq. ( |4"TT i one 
finds that 



or at least they are not much larger than their corresponding 
values for smaller absolute values of Thus, in the limit 
Ltip //i< 1 one finds for the leading behavior of the critical 
Casimir force 



JDer) _ * 
J± ~ h L ti P 



(d-i) 7+,+ 



/ 



-L>± 



1/v 



(43) 



That is, at short distances the normal critical Casimir force 
obeys an algebraic behavior L^f ^ with a power law dis- 
tinct from that for the force between two planar walls. For 
L,i p /t,± — > (or T — > T c , respectively) one finds 



JDer) 

J± 



J-LLT 
h "P 



(d-l) 



(44) 



Beyond MFT, the general trends of the scaling function for 
the critical Casimir force between planar walls are similar to 
those within MFT i 23 i 24 Therefore, we expect Eq. d43l to be 
a reliable approximation in the limit Lgp/h — > also beyond 
MFT. 

The full numerically obtained data for the normal critical 
Casimir force are in qualitative agreement with the Derjaguin- 
like approximation (see Fig. [T3T l, However, this approxima- 
tion [Eq. j43l l does not capture the variation of the amplitude 
of the force as a function of 7, which we find from the full 
numerical results [Fig. [15). 

Finally, the analysis reveals that for the normal critical 
Casimir force between a corrugated substrate and a flat sub- 
strate there is a crossover from a distant effective planar wall 



(sf) 



L ff for the case T = T c ) to a nearby regime 



"eff 



regime (/ ± 

■/tip 

which this crossover occurs is varying with the wedge an- 
gles 7: for small 7, this crossover occurs only at very small 
values of L t i p /h <C 1, whereas for large 7 even for separa 
tions L, 

nearby behavior [Fig. [TBI . 



L^T *0 [see Fig. [T3J. We find that the distance at 



tip/h ~ 1 the critical Casimir force varies according its 



C. Lateral critical Casimir force for two identically structured 
substrates 



ff er) = \ 2^rJ?S((W&=) 1/v 



up 



dL 



where f^(y±) = ^fl'± (y±). Within MFT the absolute 
values of ft'±(y±) and its derivatives with respect to y± do 
not diverge and are generally decreasing with increasing \y±\, 



If two geometrically structured surfaces oppose each other, 
in addition to the normal critical Casimir forces there are also 
lateral critical Casimir forces acting on the substrates. Lat- 
eral forces are a particularly sensitive probe of the influence 
of roughness because there is no underlying planar wall con- 
tribution. We restrict our study to the case in which both sub- 
strates are identically structured but laterally shifted with re- 
spect to each other as shown in Fig. [16] The normal critical 
Casimir force for this geometry depends on the lateral shift D, 
but its distant behavior is similar to the effective planar wall 
limit discussed in Subsec. HlIB| [Eqs. ( 139b and (l40lil. but here 
with an effective planar wall distance L e g = L — 1h e «. 

Therefore, in this subsection we focus our attention on lat- 
eral critical Casimir forces. The lateral force [Eq. (l33l l for 
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FIG. 15: Normal critical Casimir force f { l f] (L,h,y,t) for the (sf) 
geometry (see Fig. l 12t for T = T c , d = 4, and for various wedge an- 

(s f) 

gles y as a function of L t j p /h = (L/h) — 1. is rescaled by the 

factor [/z^/A+.-i-] , such that the plotted quantity is dimensionless and 
positive. All data corresponding to y = 90° are multiplied by 10 , 
the ones corresponding to y = 37° are multiplied by 10~ 2 . The gray 
lines correspond to the effective planar wall behavior [Eq. J40H with 
a decay <x L^jt, which is reached for L r i p /h ^> 1, i.e., at far distances; 

L e ff/h = L t i p /h+ 1 — h e ff(y). At small distances the actual data 
exhibit a crossover to a nearby regime with a variation oc L 



-(d-l) 
tip 



for L t ip/h <C 1. (For y= 37° this crossover occurs only for very 
small values of L,j p /h below which the numerical analysis is hardly 
feasible.) The Derjaguin-like approximation for the nearby regime 
[Eq. J44H . which does not account for different values of y, cor- 
responds to a force behavior /|° er ' x |V/A + . + ] = (Lti p /h)~( d ~ l \ 
This behavior is exhibited qualitatively by the full numerical data. 
However, the amplitude depends on y, which is described by in- 
troducing an additional factor S(y). From fits to the full numer- 
ical data we find for these factors the values S(y = 139°) = 1.0, 
S(y= 90°) = 1.4, and S{y= 37°) = 3.5; these correspond to the solid 
lines in the figure. For large wedge angles the (local) Derjaguin-like 
approximation agrees also quantitatively with the full numerical data, 
whereas for small values of y the non-local character of the critical 
Casimir force is stronger. 



this (ss) geometry exhibits the scaling form 

fW(D,L,h,y,t) = ^f{Z\D ± ,y ± ,h±,Y) (45) 



L d 

«->0 d — 1 



L d 



A +!+ A^,£, 7 ), 



where f^ s ± and Aj| are the scaling functions and the gen- 
eralized Casimir amplitude for the lateral force, D± = D/£± 
is the lateral shift scaled by the correlation length, and y± — 
t{L/^) l l v is the scaled distance variable. For T — T c the 

dimensionless quantity /,| x h d depends only on L/h, D/h, 
and y. Due to symmetry reasons it is sufficient to consider the 
lateral shift range < 8 < 1 /2, where 8 — D/l pro j. 

As will be shown below, the lateral force decays exponen- 
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FIG. 16: Confinement consisting of two identically geometrically 
.structured substrates (denoted as (ss) geometry). Both structures are 
periodic in x direction, but they are shifted by D with respect to each 
other. 



tially (i.e., stronger than L d ) with the wall distance L. This 
implies that 



'(ss) 



L>h,L, 



0. 



(46) 



The exponential decay of the lateral force as a function of the 
wall separation L for values of L that are much larger than the 
corrugation height h is expected to be present for any geom- 
etry of the corrugation. However, as a function of the lateral 
coordinate D, the behavior of the lateral force in the regime 
where the two plates approach each other depends strongly on 
the details of the geometry. One can even think of geometries 
for which the lateral force as a function of D decays alge- 
braically for D <C Ipmj- Consider the 'extreme' corrugation 
geometry of peaks appearing periodically, like a comb geom- 
etry with all teeth but every fifth tooth removed. If two such 
structured substrates approach each other such that L < 2h, 
i.e., Lfip < 0, and < D <§C l pm j, the main contribution to the 
lateral force will be due to the critical Casimir force acting 
on two vertical walls on a horizontal support. This force is 
the same as the 'normal' critical Casimir force between two 
planar walls, and it decays algebraically as a function of the 
lateral distance as long as D <C % . For D ~ lpmj/2 the influ- 
ence of the neighboring corrugation peak becomes strong and 
the behavior of the lateral force is different from the planar 
wall behavior. This consideration points out that the lateral 
force can be viewed to be similar to a normal force but acting 
on tilted surfaces. 



1. Scaling functions within mean field theory 

We calculate numerically the scaling functions in Eq. d45b 
within MFT (d = 4). We find that for a fixed distance L the 
free energy of the critical fluid in the confined (ss) geometry 
takes its minimal value at 8 = 0, i.e., the configuration of op- 
posing tips of the corrugations is the preferred one. The con- 
figuration 8 = 1/2, where ridges are opposing wedges, cor- 
responds to the maximal value of the free energy. Therefore, 
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FIG. 17: The generalized Casimir amplitude A| of the lateral criti- 
cal Casimir force [Eq. J45H for the (ss) geometry (see Fig.ll6t within 
MFT (d = 4) as a function of the lateral shift 5 =D/l pm j. The ampli- 
tudes are multiplied by \Ltip /-£-] m order to facilitate a quantitative 
comparison with the normal critical Casimir force between two pla- 
nar walls at distance L,i p . This shows that the lateral critical Casimir 
force is about a factor of 100 smaller than its normal counterpart. 
Here the wedge angle is y = 139° which corresponds to l pm j = 5.33ft. 
For L t ip < lp m j/2 the shape of the scaling function for the lateral 
force is asymmetric around 8 = 0.25, whereas for larger distances 
the shape becomes sinusoidal. 



the lateral critical Casimir force is negative (i.e., pointing in 
negative x direction) in the range < 8 < 1/2 and zero for 
8 = 0, 1/2. The lateral force is symmetric around the point 

(8 = 1/2, /jf = 0) and positive in the range 1/2 < 8 < 1. 
For the (+,+) case discussed here the 'in-phase' configura- 
tion 8 = is always the stable one. However, if the sub- 
strates are reconfigured with opposing surface fields ((+,—), 
see, e.g., Refs.|9, 10), critical Casimir forces change from be- 
ing attractive to being repulsive, and we expect the preferred 
configuration of the two substrates in the (ss) geometry to be 
the one in which the crest of an edge of one substrate opposes 
the center of a wedge of the other substrate, i.e., an 'out-of- 
phase' configuration with 5 = 1/2, in contrast to the (+,+) 
case. 

The shape of the lateral critical Casimir force as a function 
of 8 changes upon varying the distance between the walls as 
shown in Fig.[l7] Whereas for Ln p < lp m j/2 the shape of the 
generalized Casimir amplitude is asymmetric around 5 = 1/4, 
it becomes sinusoidal for L > l pro j. This asymmetry for small 
wall distances reflects the details of the geometry of the corru- 
gation, whereas the symmetric sinusoidal shape depends only 
on the major parameters of the corrugation, namely its period- 
icity and its height. 

We define the maximum of the lateral critical Casimir force 
f,i [L.h, y.t) to be its maximal absolute value in the range 

J \\,max\ tin) b 

< 8 < 1/2. The corresponding maxima of the univer- 
sal scaling functions of the lateral force are f^ s ^^(y±,h±,Y) 

and A?, (Llh,y). Typical values of A„ range from to 

||,meue v ' 1iJ JV \\,max & 

i^(100), depending on the distance between the two walls (see 
Fig- US a)). This means that the lateral critical Casimir force 



can be much stronger than the normal critical Casimir force 
between two planar walls at distance L. However, L is not 
necessarily the relevant distance to compare with (see Fig.[T6li. 
Compared with the magnitude of the normal critical Casimir 
force between two planar walls at distance L e g, the maximum 
of the lateral force can be of similar strength. The magnitude 
of the lateral force is of the order of i^(O.Ol) compared to the 
normal critical Casimir force between two planar plates at dis- 
tance L t i p as long as the corrugated walls are not too far away 
from each other (see Fig. [T8lb)). For L t j p ^> l pm j the lateral 
Casimir force vanishes relative to the normal force. Concern- 
ing the behavior for different corrugation patterns we find that 
for La p ;$ h the maxima of the lateral critical Casimir force are 
larger for small values of y (Fig. [18]T>)). On the other hand, 
the lateral force for small values of y vanishes faster upon in- 
creasing the wall distance. 

From Fig. [T8l b) we infer that the nearby behavior (for 
Ltip/h — ► 0) of the lateral critical Casimir force maximum is 

dominated by a behavior L tjp , similar to the normal 
force [Eq. 1441 1. However, there is a pronounced difference to 
the normal force concerning the distant behavior. Whereas the 
normal force approaches an effective planar wall limit with an 
algebraic decay L~j at T = T c [Fig.[l4), the lateral force de- 
cays exponentially for T ^ T c as well as for T = T c [Fig. [T9l . 
The exponential decay of the lateral critical Casimir force can 
be explained in terms of a superposition approximation for 
large distances between the walls. If one uses the approxima- 
tive forms of the order parameter profile at a single corrugated 
substrate obtained in Subsec. Ill C 21 and calculates the stress 
tensor components one arrives at (see AppendixlBt 

/ | j'"' ) (D,L>/ 1 ,/ ; ,y,f =0)~ 

- K \^A e -^ An (^D\ (4?) 

h \ 'l>roj J 

where K is dimensionless, positive, and depends only on 
y. This form captures the sinusoidal distant behavior °= 
— sin(27r5) of the lateral critical Casimir force for fixed L as 
well as the free energetic preference of the opposing tips con- 
figuration. Equation ( |47| | implies the exponential decay of the 
maximum of the lateral force: 

*^~p{-** ot (£)£}. (48) 

Although the prefactor K in the appropriate expression in 
Eq. Wl\ cannot be given analytically, the decay constants 
found using this superposition approximation agree well with 
the numerical data for L/h 3> 1 (see Fig. [l9l . The exponen- 
tial decay constant in Eq. (PPTT i corresponds to the 'reciprocal 
lattice vector' of the periodic array of wedges and ridges. It 
is worth mentioning that the lateral electrodynamic Casimir 
force due to corrugations of the walls decays exponentially 
for large wall distances with the reciprocal corrugation lattice 
vector as decay constant. 5 ^^- 

As a function of temperature the maximum strength of the 
lateral critical Casimir force is found numerically to behave 
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FIG. 18: (a) Amplitude Am n J m . (i.e., maximal absolute value with re- 
spect to the lateral shift D) of the lateral critical Casimir force scaling 
function at criticality for the (ss) geometry (see Fig.l 16b within MFT 
(d = 4) as a function of L/h for different wedge angles y. The values 

~(ss) 

of A„ indicate the strength of the lateral force compared to the 

\\.max ° r 

normal critical Casimir force between two planar walls at distance 
L (see Eq. ( 145b and Fig. 116b: their magnitudes may be similar, (b) 
Generalized Casimir amplitude A^j^ multiplied by [L,i p /L] d as a 
function of Ltip/h =L/h — 2. This product is the ratio of the maximal 
lateral force at distance L and the normal force between planar walls 
at distance L t j p (see Eq. d45t). The fact that these curves vanish lin- 

~(ss) 

early for L,j p — » (see the shaded straight line) indicates that Am 

L(jp for L,j p —¥ 0, analogous to the normal force [Eq. J44H . This 

implies that in (a) j&^L/h -> 2) °= (L/h - 2)"('^ 1 ). For small 
values of L tlp /h the lateral critical Casimir force is stronger for cor- 
rugations with a short corrugation wavelength l pro j, i- e -, at small an- 
gles y. On the other hand, the amplitudes of the lateral force for 
these short wavelength corrugations decay faster for L,i p /h 2> 1 . The 
maximum of the lateral critical Casimir force at distance L is of the 
order of ^(0.01) compared to the normal critical Casimir force be- 
tween two planar walls at distance L t i p . The curves both in (a) and 
(b) decay exponentially for L — » °° (see Fig.l 19b. 
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FIG. 19: Maximal absolute value of the lateral critical Casimir force 
[Eq. | |45H for the (ss) geometry [Fig. 1161 at T = T c rescaled by 
[/r/IA+,4-1] such that the plotted quantity is dimensionless. The lat- 
eral force decays exponentially with L/h for all wedge angles y, but 
with different characteristic decay constants. The exponential de- 
cay for L ^> h, l pro j agrees reasonably well with the predictions from 
the superposition approximation [Eq. {48}] (dotted lines). From fits 
to the curves we find for the dimensionless prefactor K in Eq. J48b 
the values K(y = 139°) =2.8 x lfr 2 , K(y =90°) = 1.5 x 10 2 , and 
K(y=3T) =7.2x 10 9 . 



and the angle y, i.e. 



L * =W i£ >y Up . 



(50) 



The numerical data for the examples presented in Fig. I20lim- 
ply that the dependence of W on Ln p /h and y cannot be further 
reduced to a dependence on L t i p /l pm j only. We find typical 
values of W to be in the range 0.8 to 1.5 for the examples 
shown in Fig. [20] However, it is rather interesting that the 
temperature dependence for the lateral critical Casimir force 
between corrugated substrates can be related to the normal 
critical Casimir force scaling function for flat substrates in a 
simple way. We interpret this aspect as a consequence of the 
lateral force to be basically a sort of 'normal force' between 
several tilted walls, generated by the same physical mecha- 
nisms. 



IV. SUMMARY 



approximately as 



(§>r) fl'M) 



(49) 



where f^(y±) = f^(y±,h± =Q,y=n) (see Eq. ©) is 
the planar wall normal force scaling function which is eval- 
uated at y± = y* ± = t (L*/^ ± ) 1 / v . This is shown in Fig. [20l 
From the numerical data we infer that the characteristic dis- 
tance L* is equal to the tip-tip distance L,j p apart from a dimen- 
sionless scaling factor W which depends on the ratio L t i p /h 



We have investigated the universal properties of strong criti- 
cal adsorption of fluids on geometrically structured substrates 
which are modeled by an array of wedges and ridges char- 
acterized by the corrugation height h and the wedge angle y 
[Fig-Hi- These universal properties correspond to the physical 
quantities of a system whenever all relevant length scales are 
larger than their characteristic molecular length scales. More- 
over, we have studied the singular contributions to the nor- 
mal and lateral critical Casimir forces acting on geometrically 
structured substrates which confine a fluid close to criticality. 
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FIG. 20: Maximum (with respect to D) of the scaling function of 
the lateral critical Casimir force for the (ss) geometry (see Fig. 1 16t 
in units of its value at T = T c [Eqs. d45t and J49H as a function 
of y* ± = f(L7<^ ± ) 1 / v (note that within MFT y* + = t(L*) 2 , and 
y* = 2t(L*) 2 ). L* is related to the tip-tip distance L llp [Fig. [16) as 
L* = W(L tip /h, y) L tip [Eq. g0>], where IV is a dimensionless scaling 
factor to be determined (see below and the main text). The gray curve 
corresponds to the scaling function for the normal force between two 
planar walls at distance L* . The data points correspond to the sub- 
strate configurations y = 90° and (i) L t i p /h =1/2, (ii) L t i p /h = 1, (iii) 
Lti p /h = 2, (iv) L,j p /h = 3, and to the configurations (v) y = 37° and 
Lti p /h = 1, as well as (vi) y = 139° and L,j p /h = 4/3. The numeri- 
cally determined scaling factors W for these cases are approximately 
(i) 1.42, (ii) 1.17, (iii) 0.94, (iv) 0.86, (v) 0.81, and (vi) 1.17 with an 
error of ca. ±0.05 for all cases. With these rescaling factors all data 
fall onto the gray curve. 



A. Critical adsorption 

Concerning the study of critical adsorption at a single cor- 
rugated substrate, in Sec. [II] we have obtained the following 
main results based on general scaling arguments and explicit 
mean field calculation for scaling functions. 

A. 1 . Close to the bulk critical point at T = T c the order pa- 
rameter profile can be described in terms of universal scaling 
functions depending on length variables scaled by the corre- 
lation length £± above (+) or below (— ) T c and the wedge 
angle y [Eq. (Q}]. Within the scaling limit the order parameter 
profile diverges close to the walls according to a power law 
[Eq. (fl2l ll, For T ^ T c it decays exponentially into the bulk 
[Eq. (3)]. At criticality the order parameter profile reduces to 
a pure power law in the direction z perpendicular to the mean 
substrate surface, multiplied by a universal amplitude func- 
tion [Eq. (O] which depends on x/h, z/h, and y, where x is 
the direction of periodicity of the corrugation. 

A. 2. The distant behavior far from the wall is governed 
by an effective planar wall behavior, i.e., the order parameter 
profile resembles the one generated by a planar wall located 
at the effective height z = h e ff [Eq. ©J with h e g — h h e ff{y) 
[Eq.©]. 

A. 3. We have calculated numerically the scaling func- 



tions within mean field theory for T < T c , T = T c , and T > T c 
[Figs.[2][3][4][8] and[7J. Approximately, the effective height is 
given by h e ff(Y) — 1 — y/2% [inset of Fig. [7J. 

A.4. In order to compare the proliferation of the sur- 
face structure into the bulk for the cases T < T c , T = T c , and 
T > T c , we have introduced the width As(s) of the undulation 
of the order parameter contour lines with mean position z = hs 
[Fig. 0. Within mean field theory we have calculated As(s) 
using a superposition approximation serving as a distant wall 
approximation. We have found that the surface structure im- 
printed onto the liquid decays exponentially into the bulk for 
all temperatures, including T c [Eq. (fT5l) l, This is confirmed by 
the numerical data. The analytically determined correspond- 
ing exponential decay constants [Eqs. ( [T6l > and ( [171 )1 compare 
well with the numerical data [Fig.[6]|. Interestingly, the lateral 
structure of the order parameter decays into the bulk fastest at 
T = T c , in contrast to the order parameter profile itself. 

A. 5. We have introduced the excess adsorption as the inte- 
gral of the difference of the order parameter and its bulk value 
over a unit cell of the periodic structure of the system. For a 
binary liquid mixture this is a measure of the total enrichment 
of a species at the substrate. Close to T c the singular part of 
the excess adsorption r exhibits scaling [Eq. d24i il with the 
corresponding universal scaling function F± depending only 
on h l%± and y. In order to compare the excess adsorption at 
a corrugated substrate with the one at a flat substrate, we have 
defined a suitably defined reduced relative excess adsorption 
[Eq. daTJ]. 

A.6. We have calculated the universal scaling functions for 
the excess adsorption within mean field theory and have com- 
pared them with the one for planar walls [Figs.l9land [T0l . We 
have found that the amount of adsorbed matter is less than 
the one at a planar wall of the same actual substrate surface. 
This effect is increasing for decreasing angles y. Compared to 
a planar substrate with the same projected area the excess ad- 
sorption on a corrugated substrate is larger. For y — > the ratio 
of these latter two adsorptions even diverges, and this effect is 
enhanced upon departing from criticality. We have defined a 
suitable ratio of the excess adsorptions above and below T c 
[Eq. (l30b l and have found that it may be smaller than the cor- 
responding planar wall ratio, depending on the wedge angle 
y. For y — > the difference between these two ratios is of the 
order of 10% [Fig.fTTI. 



B. Critical Casimir forces 

In Sec. [HI] we have investigated critical Casimir forces be- 
tween geometrically structured substrates by considering two 
arrays of wedges and ridges separated by a distance L in z 
direction and periodic in x direction [Figs. [121 and [T61 . We 
have focused on identical chemical boundary conditions on 
both substrate surfaces. In the following our main findings 
are summarized. 

B.l. Close to criticality the normal and the lateral critical 
Casimir forces obey scaling behaviors [Eqs. (|38l l and d45l)l. In 
order to determine the corresponding universal scaling func- 
tions for T < T c , T = T c , and T > T c , we have applied the 
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stress tensor method. It turns out that the improvement term to 
the stress tensor is vanishing for periodic geometries. The or- 
der parameter profiles and the universal scaling functions for 
the critical Casimir forces have been calculated numerically 
within mean field theory for the following two basic configu- 
rations. 



1. Normal critical Casimir force 

B.l.l. First, we have considered the configuration of 
a geometrically structured substrate opposing a planar wall 
[Fig- El- In this case there is a normal critical Casimir force 
acting on the substrates. We have found a crossover between 
two limiting regimes. 

B.1.2. The distant wall regime L/h ^> 1 is governed by an 
effective planar wall behavior, for which the scaling function 
of the normal force attains the form of the one between two 
planar walls separated by the distance L — h e g, where h e jj turns 
out to be the same as for a single corrugated substrate [Figs. [13] 
and[H. 

B.1.3. On the other hand, the nearby regime exhibits an 
algebraic behavior (L t i p )~( d ~ l \ where L t j p = L — h is the 
tip-wall distance [Fig. [15]. This behavior is different from the 
one between two planar walls, which would be H~ d , if H 
is the wall distance. Using a Derjaguin-like approximation, 
we have explained this behavior qualitatively [Eqs. ( f4Tb and 
<g4)]. 



2. Lateral critical Casimir force 

B.2.1. Second, we have considered the case in which two 
identically corrugated substrates are opposing each other, but 
possibly shifted laterally with respect to each other by a dis- 
tance D along the x direction [Fig. [16). The universal scal- 
ing function of the resulting lateral critical Casimir force act- 
ing on the substrates depends on D/£,±, L/£,±, h/£,±, and y 
for T 7^ T c , and on the scaled variables D/h, L/h, and y for 
T = T c , respectively. 

B.2.2. We have found that the configuration of opposing 
tips is the preferred one and the configuration of a wedge op- 
posing a tip corresponds to an unstable configuration. 

B.2.3. The shape of the scaling function of the lateral force 
in between those two configurations depends on the relative 
shift D between the walls [Fig. [IT)- For wall distances which 
are small compared to the projected width of a wedge or the 
height of the wedges the lateral force reflects the detailed ge- 
ometry of the system. On the other hand, for distances larger 
than h and the projected width of a wedge the shape becomes 
sinusoidal [Fig. [171 . 

B.2.4. The amplitude of the lateral critical Casimir force 
depends on L but there is also a strong dependence on the 
wedge angle y [Fig. [18). Within a suitable comparison scheme 
the amplitude of the lateral critical Casimir force turns out to 
be of the order of 1% of the one for normal critical Casimir 
forces between two planar walls at distance L — 2h. There is 
evidence that in the nearby regime the lateral force amplitude 



scales as (L — 2/i)~^~ 1 ) [Fig. [T8l b)l. The distant behavior 
is characterized by an exponential decay of the lateral force 
amplitude into the bulk [Fig. [19). 

B.2.5. Using a superposition approximation, we have been 
able to explain the sinusoidal shape of the lateral force, the 
free energy behavior, and the exponential decay of the force 
amplitude for large distances [Eq. (07)] ■ The approximatively 
calculated decay constants compare well with the full numer- 
ical data [Fig. [19). 

B.2.6. We have found that the temperature dependence of 
the maximal strength of the lateral critical Casimir force can 
be expressed in terms of the temperature dependence of the 
normal critical Casimir force between planar walls [Eq. ( 149) 
and Fig. [20). This can be interpreted in the sense that the lat- 
eral critical Casimir force between corrugated substrates can 
be viewed to be a normal critical Casimir force between tilted 
walls, and thus being generated by the same physical mecha- 
nisms as the normal critical Casimir force. 



3. Comparison with the quantum electrodynamic Casimir effect 

Surface corrugations lead to pronounced ef- 
fects on the quantum electrodynamic Casimir force, 

t00 4 9.50.51.52.53.54.55.56.57.58.59.60.61 

B.3.1. Walls, which are geometrically structured, lead 
to a behavior of the normal quantum electrodynamic 
Casimir effect, which is different from the one for planar 
walls i 49 ' 50 ' 51 ' 52 ! 53 Generally, the normal quantum electrody- 
namic Casimir force is enhanced compared with the one be- 
tween two planar walls separated by the same mean distance, 
and weakened compared with the one between two planar 
walls at the same minimal, i.e., tip-to-tip distance. This is 
similar to the behavior of the normal critical Casimir force 
[Fig- El a)], which we also found to be enhanced compared 
with the one acting between planar walls separated by the 
same mean distance, and weakened compared with the one 
acting on planar walls at distance L r , p [Fig. El . 

B.3.2. The effective planar wall limit of the normal crit- 
ical Casimir force [Fig. [T3l b)l with the effective corruga- 
tion height h e ff = hh e ff(y) given by the expression h e ff(y) for 
semi-infinite systems [Eq. (120) 1 is consistent with the results 
obtained for the normal quantum electrodynamic Casimir 
force between a periodically structured and a planar wall. 51 
However, most corresponding studies for the quantum elec- 
trodynamic Casimir effect consider sinusoidal or rectangu- 
lar shaped corrugation profiles, which are different from the 
wedge-like shapes which we have discussed upon extending 
previous studies of this type of geometry i 46 ' 47 Therefore, a 
comparison of the nearby regimes of the two effects is diffi- 
cult. 

B.3.3. The dependence of the amplitude of the lateral criti- 
cal Casimir force for the (ss) geometry on the wall distance 
[Fig. E), with a geometry dependent form in the nearby 
regime and a sinusoidal form in the distant regime, is simi- 
lar to the one for the lateral quantum electrodynamic Casimir 
force i 5 - 5 - 
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B.3.4. Both the lateral critical Casimir force for geometri- 
cally structured substrates [Fig. [TTJ and the lateral quantum 
electrodynamic Casimir force are directed such that the pre- 
ferred configuration of the two corrugated substrates is the one 
of opposing tips i 51 ' 54 ' 55,56 However, critical Casimir forces 
change from being attractive to being repulsive, if the sub- 
strates are endowed with opposing surface fields ((+,—), see, 
e.g., Refs. 10.19.22.24) instead of the (+,+) case discussed 
here. For corrugated substrates corresponding to the (+,— ) 
case we expect the preferred configuration of the two sub- 
strates in the (ss) geometry to be the one in which the crest 
of an edge of one substrate opposes the center of a wedge of 
the other substrate, in contrast to the (+, +) case. 

B.3.5. The exponential decay of the lateral critical Casimir 
force amplitude for large distances (see Eq. (|48) ) is consis- 
tent with the behavior of the lateral quantum electrodynamic 
Casimir force , 51 ' 56 too. Both decay with the 'reciprocal lattice 
vector' corresponding to the periodic lateral structure as the 
characteristic inverse decay length. 

B.3.6. The nearby regimes of the lateral critical Casimir 
force and the lateral quantum electrodynamic Casimir force 
are difficult to compare, because most of the corresponding 
studies of the lateral quantum electrodynamic force deal with 
sinusoidally shaped corrugations, whereas the present study 
of the lateral critical Casimir force corresponds to wedge-like 
shapes. 

B.3.7. Actuation of nano-mechanical devices by the quan- 
tum electrodynamic Casimir effect may be possible using geo- 
metrically structured objects i 58 ' 59 ' 60 ' 61 We expect that similar 
applications may be achieved using critical Casimir forces. 
However, the fact that the critical Casimir force can be 
changed from being attractive to being repulsive upon chang- 
ing the chemical boundary conditions opens up additional 
possibilities. 



C. Discussion and outlook 

C.l. In addition to the critical Casimir forces due to fluctu- 
ations of the order parameter, there are the omnipresent van 
der Waals dispersion forces acting as effective background 
forces on the substrates confining the fluids. The total force 
is approximatel y 77 ' 78 the sum of these background forces 
and the critical Casimir force. However, it has been shown 
theoretically 28 and experimentally (see, e.g., Ref.l22l) that for 
classical fluids (corresponding to the presence of symmetry 
breaking surface fields) the normal critical Casimir forces for 
T — > T c are much larger than these background forces and 
are directly accessible. The same is expected to hold for the 
normal critical Casimir force between corrugated substrates. 
We have derived an estimate for the strength of the lateral 
critical Casimir force compared with the strength of the lat- 
eral dispersion background forces for such systems (see Ap- 
pendix0. Accordingly, for typical physical parameters [Sub- 
appendix IC 21 and upon approaching criticality the lateral crit- 
ical Casimir force is much larger than the corresponding lat- 
eral van der Waals force. 

C.2. Lateral critical Casimir forces also occur if simple 



fluids or liquid crystals in their nematic phase— are con- 
fined by chemically structured substrates. Also in these sys- 
tems the lateral force is sinusoidally shaped for large substrate 
distances and depends on the details of the pattern structure 
for small substrate distances. Similar to the lateral critical 
Casimir force for the (ss) geometry (see Fig. [19} also in these 
systems the lateral forces decay exponentially for increasing 
wall distance with the 'reciprocal lattice vector' as character- 
istic length scale. 44 The strength of the lateral critical Casimir 
force due to the confinement of critical fluids by two sub- 
strates patterned with chemical stripes is of the order of the 
normal critical Casimir force between these substrates when- 
ever the wall distance is small. 37 This ratio is larger than the 
corresponding one of the lateral critical Casimir force for the 
(ss) geometry compared to the normal critical Casimir force 
between two plates at distance L tlp (see Figs. [T6land [T8l b)). 
However, if one compares the lateral critical Casimir force 
with the normal one between two planar walls of distance L 
(see Figs. [T6land [T87 a)) or L e ff, one finds that its strength can 
be of similar order (for L e jf) or even larger (for L) compared 
with the corresponding planar wall normal force. Therefore, 
lateral critical Casimir forces due to corrugations may be com- 
parable in magnitude with lateral critical Casimir forces due 
to chemical structures. 

C.3. An experimental realization for the measurement of 
normal and lateral critical Casimir forces acting on geomet- 
rically structured substrates would provide access to the ge- 
ometrical structure of the corrugation of buried solid-liquid 
interfaces. Recently, the direct measurement of the critical 
Casimir force between a colloid and a planar wall confining 
a binary liquid mixture has been reported. 22 A similar experi- 
ment, in which the flat substrate is replaced by a corrugated 
one, with a typical corrugation size much smaller than the 
extension of the colloid, could provide access to the normal 
critical Casimir force in the (sf) geometry (see Figs. l2TT a) 
and [12). The lateral quantum electrodynamic Casimir force 
has been measured between a sinusoidally corrugated metal 
sphere and a corrugated planar wall. 54 An experiment involv- 
ing a geometrically structured colloid could provide direct ac- 
cess to the lateral critical Casimir force for the (ss) geometry 
(see Figs. l2TT b) and [Toll. Such setups can play a role for mov- 
ing ratchet-like parts in nano-machines. Using a smooth col- 
loid which is small compared with the corrugation, one may, 
for example, measure the individual local contributions to nor- 
mal and lateral critical Casimir force [Fig.l2TTc)l. 
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FIG. 21: Possible experimental realizations for the measurement 
of critical Casimir forces for geometrically structured confinements 
(schematic). Commonly, in experiments colloids or spheres are used 
instead of two perfectly aligned planar walls (see, e.g., Refs. l22ll54ri . 
Here the colloids are considered to be long rods, (a) A 'large' colloid 
opposing a geometrically structured substrate could provide access 
to the normal critical Casimir force for the (sf) geometry fFig. 1121 . 
(b) A colloid with imprinted lateral sinusoidal corrugations opposes a 
geometrically structured substrate (a similar setup has been realized 
for the quantum electrodynamic Casimir force). 54 In such a geom- 
etry, lateral critical Casimir forces should arise similar to the ones 
arising in the (ss) geometry fFig. 1 1 61 . (c) A 'small' colloid which 
is comparable in size with the corrugation height and the lateral ex- 
tension of a single wedge could provide access to the normal or the 
lateral critical Casimir force acting on a single 'tip' of the upper sub- 
strate in Fig. 1161 



APPENDIX A: UNDULATION OF THE ORDER 
PARAMETER CONTOUR LINES 

The behavior of the undulation As(s) of the contour lines of 
the order parameter profile close to a single corrugated sub- 
strate (see Subsec. HlC21 i is studied within MFT and a distant 
wall approximation scheme by introducing a superposition ap- 
proximation for the order parameter [Eqs. (\3[ and (I14H , 



1. Distant behavior of the order parameter profile for T — T c 



and the fact that the underlying effective planar wall contribu- 
tion is a solution of the corresponding Euler-Lagrange equa- 
tion Eq. ( fTTT i, Eq. ( IA1 b turns into 



<\x) 



mf{z) 



- 3 m, 



(«?>«)' 



(A4) 



Equation ( IA4t holds for x and z independently and therefore 
both sides of Eq. (lA4b are equal to the same constant k 2 so 
that for the part of Eq. iA4i depending on x one has 



dx 2 ' 



nf{x) = -tfmf{x) 



(A5) 



A solution of Eq. (IA5t . which is compatible with the periodic 
boundary conditions of the corrugated system [Fig.Q~|, is 



(x) 

m Q (x) = —M x cos(kx) 7 



(A6) 



where k = kj^- with k a natural number; M x is a positive am- 

Iproj 

plitude. We assume that far from the substrate only the main 
contribution due to the corrugation is significant, i.e., k = 1 
and k = 2n /l pro j. For the z-dependent part of Eq. (IA4l i one has 



< J (z). 



(A7) 



Far from the substrate (i.e., z lproj) the underlying effective 
planar wall contribution Eq. (I A3t fulfils 
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1 l proj 



3 V/, 



(A8) 



Thus, together with k = 2%/l pro j from above, Eq. (IA71 > can be 
written as 



Within MFT the order parameter profile 
m(x,z,h, y,t = 0) = mo(x,z) at criticality solves the dif- 
ferential equation [Eq. ( fTTT ll 



V 2 mo(x,z) = (m (x,z)) 3 



(Al) 



Far from the substrate, i.e., for z/h ^ 1, the order parame- 
ter profile contribution mk\x,z) due to the corrugation and 

the underlying planar wall contribution (z) [Eq. (fT3l l are 
related as 



(A2) 



Thus, for the cubic power of mo one can neglect quadratic 

or higher order terms of m^. Within MFT the underlying 
effective planar wall contribution [Eqs. © and ([8]l] is 



m o (z) = c+{z- 



where c + is the universal surface amplitude for the planar wall 
system [Eq. ©J with c+ = V2 within MFT. 15 Together with 

Eq. ( TBi i (i.e., the product ansatz m^^Xjz) = (x) (z)) 



with the physically relevant solution 



o W W = 



, 27t 

M z exp<j -- — z 

lproj 



(A9) 



(A10) 



and M z as a positive amplitude. Accordingly, at criticality the 
contribution of the corrugation to the order parameter decays 
exponentially into the bulk, in contrast to the total order pa- 
rameter, which decays algebraically. Upon inserting the ef- 
fective planar wall solution, the final approximative expres- 
sion for the order parameter reads 



m (x,z) 



My, COS 



\ l pro] J 



(All) 



(A3) where M xz — M X M Z . M xz depends only on the corrugation 



height (0= 1 /h) and on the wedge angle 7 because for T = T c 
the product m x h depends only on x/h, z/h, and 7 [Eq. ©]• 

From Eq. ( IA1 II ) we calculate the undulation As(s) of the 
contour lines, which has been defined in Subsec. HlC2l as the 
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difference of the maximum and the minimum value of z/h of 
a given contour line of the order parameter. The positions of 
the maximum and the minimum z max = h(s + jAs) and Zmin — 
h(s — jAs) of a contour line correspond to the same x values 



as the crest of an edge or the apex of a wedge of the periodic 
array of wedges and ridges, i.e., x = l pr0 j/2 oix = 0. Thus, for 
any contour line far from the substrate one finds 



m o( x — Iproj/'^iZ — Zmax) — c + 



h[s + -As-heff(y) 



-M xz e 'i' m r 1 



(A12) 



and 



rno(x = 0,z=z m m) = c+ 



h[s--As-heff{y) 



-l 



■M X7 e w 



(A13) 



Along a contour line the order parameter value is constant 
which allows one to equate the right-hand sides of Eqs. ( IA12b 
and ( TABl i: 



2M xz h 

2kI ' s 



cosh ( ^— As 



As 



(A14) 

Far from the wall As(s) <C 1, and quadratic and higher order 
terms of As(s) can be neglected. With the additional assump- 
tion that h is not much larger than l pm j, one can approximate 

cosh ( j^-As ) ~ 1 and Eq. ( 1A14I ) becomes 

V 'Proj J 



„M xz h 
As{s) = 2 1 s 



2 _2*!l, 
heffir)} e w 



—{s-heff(r)) 



(A15) 



Thus, the undulation of the contour lines of the order parame- 
ter decays exponentially into the bulk with the decay constant 



°o{Y) = — =7TCOt(7/2). 

'proj 



(A16) 



2. Distant behavior of the order parameter profile for T > T c 



which corresponds to the bulk value (z+ — ► °°) = 0. Far 

from the substrate P±' <C P+ , so that one can neglect terms 

of quadratic or higher order in P+\ Using the product ansatz 

according to Eq. (i.e., p[ 1) (x+,z+) = Pj* ) (*+)P+ ) (z+))> 
upon inserting Eq. ( IA1 8b into Eq. dA17b one finds the follow- 
ing differential equations: 



^\x + )=k 2 P { *\x + ) 



and 



-l-3(P? ) (z+)) 2 = k 2 



(A20) 



(A21) 



Similar to Subappendix IA li the solution of Eq. dA20l ) is 

-P+ cos(kx+) (A22) 



where P+ > and k = 2kn£,+ / l pr0 j with a natural number k. 
By keeping only the main contribution k=l, Eq. ( IA2 1 b be- 



The order parameter profile for T > T c is characterized by 
the scaling function P + (x + , z+ , h + , y) [Eq. (Q}] which we ab- 
breviate as P + (x + ,z+). The procedure follows the one of Sub- 
appendix [AJ] for T = T C . The MFT Euler-Lagrange equation 
for the scaling function reads 

(A17) 

According to Eq. ([T3l one separates P + into an effective planar 
wall contribution and a corrugation contribution: 



P + (x + ,z + )^Pl°\z + )+Pl L >(x +J z + ), 



(i) 



(A18) 



where P+ (z+) is the effective planar wall MFT solution (see 
Eqs. (|3]-[8)andFig.[D: 

pj 0) (z+) = V2/ S mh(z + -h+h eff (Y)), (A19) 



d_ 2 

dz 




P| z) (z+) 



(A23) 



where P+ (z+) is replaced by its bulk value zero because 

P+ (z+) < 1 for large distances z > h,%+- Equation ( IA23b 
is solved by 



P^(z + )=P+" 



(A24) 



where P z is a positive constant. Finally, for z+ ^> 1 the ap- 
proximate solution for P + reads 



P + (x +1 z + )^P { °\z)-P:P; 



(A25) 
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with given by Eq. ( IA19I ). For z+ ^> 1 one has 

H 0) (z+»l)-exp(-z + ). 



(A26) 



Similar to Subappendix I A II one can directly identify the 
minimum and the maximum of the contour lines with (x + = 

, Z+min = h+(s-jAs)) and (x + = ^Ipro) , Z+max = h + (s + 



I As)). Using the abbreviation 



where PT > and k = -r-^. Accordingly, the part of 

_ l pwj 

Eq. JA32b depending on z_ reads 



d 2 



j>W(z_) 2 2 



Far from the wall one can replace P_' (z_ ^> 1 ) by 1 : 



(A34) 



Q = h+\ 1 



V 4 



pro; 



(A27, §(l*V»»-£ + (*£ 

^ z V l proj 



the comparison of the order parameter at these two positions 
yields 

Ze- Qs cosh (§ Ai ) = e ~ h+s sinh (y Ai ) . ( A28 > 



(A35) 



With this replacement Eq. (IA32I ) is the same as Eq. ( IA23I ) for 
T > T c , with E,- replacing §+. Thus, with Eq. <(4j the undula- 
tion for T < T c can be approximated by 



As(s) ■ 



(A36) 



where Z is a prefactor, which is independent of s and includes where 

all possible other prefactors. Far from the substrate As <C 1, 

and one can neglect terms of quadratic or higher order of it. 

With the additional assumptions that h + is not much larger O (/?-,/) 

than unity and l pro j is not much smaller than t, + (i.e., Q is not 

much larger than 1) Eq. (IA281 > simplifies and one finds 




(A37) 



As(s) °= exp{-(Q-h + )s}, 



(A29) 



where the proportionality constant depends on h + and y. 
Thus, for T > T c the undulation decays exponentially into the 
bulk with the decay constant 



1 I —ncot(- 



(A30) 



3. Distant behavior of the order parameter for T <T C 



APPENDIX B: DISTANT BEHAVIOR OF THE LATERAL 
CRITICAL CASIMIR FORCE FOR T = T c 

In this appendix we discuss the distant wall behavior of the 
lateral critical Casimir force between two identically struc- 
tured substrates (see Subsec. IIII C\ by using a superposition 
approximation based on Eqs. ( fT3l and ( fT4t and AppendixlAl 



In Subappendix IA 1 I the corrugation contribution m Q 



to the 



order parameter too ~ m Q + m Q for T — T c close to a single 
corrugated wall [Fig.[T| has been approximated for large dis- 
tances from the wall as [Eq. (IA1 lH 



For T < T c one has to consider the scaling function 
P_(x_,z_,/i_,y) [Eq. ([I])] abbreviated as P_(x_,z_). It obeys 
the differential equation [Eq. (fTTI ll 



/>_(x_,z_) , Pl(x-,z-) 



(A31) 



Following an analogous procedure as in Subappendices IA 1 
and [A2] and using the ansatz P = pW+pWp® [ Eqs (TJJ, 
and (fl~4T >1 one finds 



dx_ 

~p¥[x~) 



P W (z_) 



The part of Eq. ( IA32I ) depending on x_ is solved by 



f W (x_) = -P~ cos(fcc_), 



(A32) 



(A33) 



- M xz cos 



2k 



(Bl) 



where M xz is a positive number which only depends on h as 
h and y. 

For the (ss) geometry [Fig. [TBI one has to deal with the 
corrugations of both walls located at z — and z — L as well 
as with the lateral shift D along the x direction. The corru- 
gation contribution of the lower and upper wall in Fig. [16] 

to the order parameter is to^^z) [Eq. ( IB 11 1 and m^\x — 
D,L — z), respectively. The underlying planar wall contri- 
bution for the (ss) geometry is the order parameter profile 

TOq ' = TOq ' + (z,L e ff) between two planar walls at distance L e g 
in the strong adsorption limit for identical chemical boundary 
conditions 10 ; m^' + (z, L e £) does not depend on the lateral posi- 
tion variable x. In order to determine the distant wall behavior 
it is convenient to approximate the total order parameter at 
T = T c in the (ss) geometry as 
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m (x,z,L > h,D,h,y) ~ ' + (z,L e ff) + m^\x,z) + m^\x - D,L — z) 



^(z,LS -M X A cos 



2k 



cos( &-(x-D) 



-(L-z) 



(B2) 



The lateral critical Casimir force is given by the integral over the stress tensor component T xz (see Eqs. d34l i and (|36ll): 



& s) (D,L,h, 7 ,t=0) = -- 
» u I, 



3! 



dx 



dmo(x,z,L,D,h,Y) 



dz 



dmo(x,z,L,D,h,Y) 



dx 



(B3) 



where zo is an arbitrary but fixed position in between the two 
substrates. For simplicity we choose zo = L/2, i.e., the center 
position in between the two substrates, where the planar wall 
contribution ' to the order parameter takes its minimal 
value, i.e., 



~{z,Leff) 



dz 



dm^ + {z 1 L ej 
dx 



0. 



(B4) 



With this the integrand in Eq. ( IB3b reduces to a product of 
the derivatives of the sum of the corrugation contributions. 
Performing the integration in Eq. ( |B3t leads to Eq. ( |47] i with 
the dimensionless amplitude 



? 3! , 
K = \% 2 — A 4 
it 



i-i 



Mjh d 
I 2 ■ ' 

pro] 



(B5) 



where d = 4 within MFT. Note that l pmj = 2/itan(y/2) and 
M xz is inversely proportional to accordingly, K is positive 
and depends only on y. 



APPENDIX C: COMPARISON OF THE LATERAL 
CRITICAL CASIMIR FORCE WITH LATERAL VAN DER 
WAALS FORCES 

1. Lateral van der Waals forces for geometrically structured 
confinements within the approximation of summing pair 
potentials 

Van der Waals or dispersion forces provide a nonsingular 
background contribution which adds to the critical Casimir 
forces. For example, the particles constituting the confining 
substrates in the (ss) geometry [Fig. Q2) do not only interact 
with the fluid, but also among each other, such that their inho- 
mogeneous distribution gives rise to a lateral force acting on 



the substrates even in the absence of the fluid in between them. 
The strength of these background forces can be determined by 
analyzing the confining system without a fluid. 

In order to estimate the contribution of the background in- 
teraction between the substrate molecules to the total force, 
we assume their corresponding pair potential to be given by 
the well known non-retarded van der Waals potential, 



«(r,r') 



1 



n 2 pp' |r- 



-'16' 



(CI) 



where r = (x,y,z) and r' = (x' ,y' ,z') are the three-dimensional 
position vectors of the two interacting molecules, the ampli- 
tude srf characterizing the strength of the interaction between 
the molecules is chosen such that it turns into the Hamaker 
constant, and p and p' are the number densities of the par- 
ticles in the two interacting bodies (here p' = p because we 
consider identical substrate materials). 

We use the pairwise summation approximation (PWS) in 
order to calculate the total van der Waals interaction potential 
between two macroscopic bodies, 

U vdw^_ f r d 3 rd*r> pp'uirj), (C2) 
Jb Jb' 

where B and Bl denote the volumes of the two interacting sub- 
strates. Note that for the geometries under consideration, the 
typical sizes of which are large on molecular scales, assum- 
ing a non-retarded van der Waals potential as well as the PWS 
amount to a severe approximation. But here using them is ap- 
propriate because we want to obtain a simple comparison be- 
tween the strengths of the lateral van der Waals force and the 
lateral critical Casimir force for the (ss) geometry [Fig. [16). 

The van der Waals interaction potential per area A = l pro j x 
l y within the xy plane corresponding to the lateral van der 
Waals force for the (ss) geometry is given by 



jjvdW ly Zmta( x . 



— s-r dx dy dz dx' dy' / dz' T) (C3) 

n 2 A J J J J J J (t x -x>Y + (y-y'Y + {z-z') 2 f 
-L mJ /2 o o — — z m ,„(x') u ) v y ) \ ) ) 
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where z m ax(x) is the z value as a function of x of the lower substrate surface in Fig 

( 2/i i| 



Zrnax (x) 



2 — — 2 ' 



^Zmax(x-^- lproj)i X ^ zb : 



(C4) 



and Zmin(x) is the corresponding z value of the upper substrate surface. For identical but shifted corrugations one has z m ,„(x) = 
L — Zmux{x — D), where D is the lateral shift between the two substrates along the x direction [Fig. Q2). In Eq. dC3b only the 
volumina of the corrugations on top of the two planar substrates have been considered. Jj vdw — U" lw is the interaction potential 
between two planar semi-infinite substrates at distance L; it depends only on z and therefore does not give rise to a lateral force. 
With P^dy{a 2 +y 2 )~ 3 = 37r/(8a 5 ), the integral in Eq. (IC3l reduces to 



vdW 



lproj/2 Zmax{x) 

dx I dz I dx 1 



dz' 



1 



((x-x') 2 + (z-z') 2 ) 5/2 ' 



(C5) 



-lpwj/2 -oo L-Zmatf-D) 

The lateral van der Waals force f^ dw per area corresponds to the derivative of /A with respect to the lateral shift D, 

vdW _ h"°jl 2 Zmax(x) 



fVi 
J 



vdW 



d u. 



dD A 



7>si 



dx dz dx' 



Z max(x D) 



lproj/2 



((x-x') 2 + (z + z max (x> -D)-L) 2 f 2 ' 



(C6) 



where z' m ax(x) = zj^Zmax(x)- Rescaling the lengths in the integrals in Eq. (lC61 l with the projected width l pm j, 

X=x/lp ro j : X —X /Iproji Z—z/lproj^ 8—Djlp r oj : L = L/lp ro j, 

and defining the dimensionless function z m ax (corresponding to Eq. ( IC4b and with h = l pm jCOt(y/2)/2 [Fig.Q~|) as 



Zmax{x 7 Y) '■ 

the lateral van der Waals force [Eq. (IC6H reads 



cot( 7 /2)|x|, 

Zmax ~F 1 ) ) 



3j^ / 1 



4<*<i 

x^ ±- 
x < =c 2' 



(C7) 



(C8) 



l7l \ Iproj 



f7 w (8,L,Y), 



(C9) 



where /it plays the role of a dimensionless scaling function of the lateral van der Waals force for the (ss) geometry: 



fi dW (5,L,Y) 



1/2 Zmax(x,y) 

J dx J dz I tlx' 

1/2 



Zmax(x'-8,Y) 



([x-x') 2 + {z + z m ax{x' -5,7)-L) 2 ) 



5/2' 



(CIO) 



With z'maxfcY) = 4zZmax{x,Y) 



r 



2. Numerical comparison 

In order to gauge the strength of the lateral critical Casimir 
force [Subsec. [ill CL we compare its amplitude with the one 
of the lateral van der Waals background force obtained within 
the PWS approximation [Subappendix IC 111 for a specifically 
chosen geometry. 

For the amplitude of the lateral critical Casimir force for 
the (ss) geometry we consider at T = T c the specific configu- 
ration 7 = 90° and L/h = 4 (i.e., L tip /h = 2) [Fig. [16). For this 
choice the corresponding lateral critical Casimir force scaling 



function amplitude has approximately the value (see the solid 
curve with square symbols in Fig. fTSl b)) 

^llj L / h = 4,7-90°) x [L tip /L] d ~ 4.1 x lO" 3 . (Cll) 

The corresponding maximum of the lateral critical Casimir 
force is (see Eq. d45l l) 



f 



d-l 

hip 



Xl/Ky) 



(C12) 



The lateral critical Casimir force /,| is given in units of 
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the characteristic thermal energy ksT c and the d— 1 dimen- 
sional cross-sectional area A within the xy plane. Accord- 
ingly, the maximum of the lateral critical Casimir force per 
area f^'^JJ-'/h,"/) at criticality for the configuration chosen 
above is given by 

flZ^/h = 4, 7 = 90°) = k B T c f^ ax . (C13) 

Although the scaling function was determined within 
mean field theory, its form gives a reasonable approximation 
for dimensions d < d uc — 4, too, if one replaces all quanti- 
ties depending on the dimension of the system by their corre- 
sponding value for the actual dimension d. In order to obtain 
an estimate of the strength of the lateral critical Casimir force 
in d = 3 we use |A+ +| = 0.44. 24 As typical length and tem- 
perature scales we use L = 200 nm (i.e., h = 50 nm, L„ ; , = 
lOOnm, and l pro j = lOOnm for the specific geometry under 
consideration) and we take T c = 300 K. With these values, the 
maximum of the lateral critical Casimir force [Eq. ( 1C13H is 
of the order of ff^ ~ 15 fN/(jUm) 2 . 



In order to estimate the strength of the contribution of the 
background forces to the total force we use the PWS approx- 
imation result for the lateral dispersion force for the (ss) ge- 
ometry [Subappendix IC II . For the lateral van der Waals force 
scaling function for the geometry defined above (Ljh = 4, 
Y = 90°, i.e., L = L/l pmj = 2) one has [Eq. dCTOH 

7^(5 = 0.25, 1 = 2, 7=90°) ~ 1.3 x 10~ 4 . (C14) 

Typical values of the Hamaker constant srf are of the order of 
IQkgT . Using l pr0 j = 100 nm as above, we find for the maxi- 
mum of the lateral van der Waals force [Eq. (IC9H for the (ss) 
geometry a typical value of l/J^J — 0.64fN/(jUm) 2 . This 
means that the background van der Waals interaction calcu- 
lated within the PWS approximation, which represents the 
background interaction, is of the order of 1% to 10% of the 
lateral critical Casimir force. Thus, near criticality the lateral 
critical Casimir forces are expected to be much stronger than 
the background forces. 
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